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Abstract 

This article reviews recent developments in statistical field theory far from equi- 
librium. It focuses on the Kardar-Parisi-Zhang equation of stochastic surface growth 
and its mathematical relatives, namely the stochastic Burgers equation in fluid me- 
chanics and directed polymers in a medium with quenched disorder. At strong 
stochastic driving ~ or at strong disorder, respectively - these systems develop 
nonperturbative scale-invariance. Presumably exact values of the scaling exponents 
follow from a self-consistent asymptotic theory. This theory is based on the concept 
of an operator product expansion formed by the local scaling fields. The key differ- 
ence to standard Lagrangian field theory is the appearance of a dangerous irrelevant 
coupling constant generating dynamical anomalies in the continuum limit. 



Contents 



1 Introduction 








1.1 Directed growth, Burgers equation, and polymers 




1.2 Overview of this article 



2 The field theory of directed strings 

77L Correlation functions and the operator product expansion 

2.2 Renormalization 

2.3 Results and discussion 



3 Directed strings in a random medium 

3.1 Replica perturbation theory 

3.2 The strong coupling regime 

3.3 A string and a linear defect 

3.4 A string and a wall 

3.5 Strings with mutual interactions 

3.6 Upper critical dimension of a single string 
JTF Discussion 



4 Directed growth 



4.1 


Field theory of the Kardar-Parisi-Zhang model 


4.2 


Dynamical perturbation theory 




4.3 


Renormalization beyond perturbation theory 





O Response functions in the strong coupling regime 



4.5 Height correlations in the strong coupling regime 

4.6 Dynamical anomaly and quantized scaling . . . 

4.7 Discussion 



Appendix 



Bibliography 



1 Introduction 



The last 25 years have seen a very fruitful application of continuum field theory to statis- 
tical physics A system in thermodynamic equilibrium tuned to a critical point "looks 
the same on all scales", that is, it shows an enhanced symmetry called scale invariance. 
Since there is no characteristic scale, the correlation functions in a critical system can only 
be power laws. Microscopically quite different systems share certain aspects of their criti- 
cal behavior which are therefore independent of the small-scale details. This phenomenon 
is called universality. Among the universal features are the exponents characterizing the 
power laws of correlation functions. 

The renormalization group has provided the first satisfactory explanation of universal- 
ity as well as calculational tools to compute, for example, critical exponents. The central 
concept is that of a flow on the "space of theories". In equilibrium systems, this flow is 
parametrized by the coupling constants of an effective Hamiltonian determining the parti- 
tion function. It can be understood as a link between the microscopic coupling constants 
defining a model and the effective parameters governing its large-distance limit. (For an 
Ising model, e.g., the former are the spin couplings on the lattice and the latter are the 
reduced temperature and the magnetic field.) The fixed points of this flow represent the 
universality classes. At a fixed point, a system can be described by a so-called renormal- 
ized continuum field theory, that is, a theory where all microscopic quantities such as the 
underlying lattice constant no longer play any role. The universal properties can often be 
calculated in a systematic perturbation expansion. 

About a decade after Wilson's seminal work on renormalization, the advent of con- 
formal field theory triggered a new interest in two-dimensional critical phenomena 0. It 
was realized that these systems have an even larger symmetry called conformal invariance. 
(A conformal transformation is any coordinate transformation that locally reduces to a 
combination of scale transformation, rotation and translation. The rescaling factor is now 
allowed to vary from place to place.) Conformal invariance severely constrains the struc- 
ture of continuum field theories in two dimensions, yielding a (partial) classification of the 
universality classes and their exact scaling properties. Under some additional conditions, 
there is only a discrete set of solutions: the possible values of the critical exponents are 
quantized. 

The subsequently developed S-matrix theory has extended the exact solvability to 
theories close to a critical point, i.e., with a large but finite correlation length. For the 
first time, it has been possible to verify the ideas of scaling beyond perturbation theory 
for a whole class of strongly interacting field theories. On the other hand, the conformal 
formalism entails a shift of focus from ^f/o 6a/ quantities (such as the Hamiltonian of a theory 
and its coupling constants) to local observables, i.e., the correlation functions. Their 
structure is encoded by the so-called operator product expansion, a rather formidable 
name for a simple concept (explained in Section 2). The operator product expansion 
is at the heart of a field theory. Conformal symmetry imposes a set of constraints on 
the operator product expansion which lead to expressions for the correlation functions. 
This theory makes no reference to a Hamiltonian. Without a renormalization group flow, 
however, the link to microscopic model parameters is lost. Hence, it has to be determined 
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Fig. 1: (a) STM snapshot of a kinetically roughened gold film. The projected area of the sample 
is 510 X 510nm^. Courtesy of Herrasti et al. js). (b) Part of a discretized KPZ surface. The 
entire surface has 400 x 400 lattice points and is shown after 1000 time steps. 

a posteriori which lattice models belong to the universality class of a given conformal field 
theory. 

Another important shift of focus has taken place in recent years. Traditionally scale 
invariance has been associated with second order phase transitions, which requires fine- 
tuning of the model parameters to a critical manifold. However, it became clear that there 
are many systems that have generic scale invariance in a region of their parameter space. 
Perhaps the simplest such systems are interfaces. For example, the surface of a crystal 
in thermal equilibrium fluctuates around an ideal symmetry plane of the crystal. The 
displacement can be described by the "height" field h{r), the two-dimensional vector r 
denoting the position in the reference plane. Above a certain temperature, the surface is 
rough, i.e., it develops mountains and valleys whose size is typically some power of the 
system size. The correlation functions of the height field become power laws as well. The 
roughness turns out to be even stronger if the crystal is growing, which drives the surface 
out of equilibrium. There is another important difference between equilibrium surfaces 
and driven surfaces. In the former case, the height pattern shows an up-down symmetry 
between mountains and valleys. Out of equilibrium, that symmetry is lost. 

Obviously, such open systems are quite common: any growth, pattern formation or 
reaction process propagating through some localized boundary or front generates a driven 
interface, often with long-ranged correlations P, ^ . An example is the kinetic roughening 
of thin metal films by vapor deposition. Fig. 1 (a) shows the ST microscope analysis of 
a gold film at room temperature from which the authors were able to extract power law 
behavior of the height correlations, indicative of a scale-invariant surface state 0. This 
state is seen to be directed (i.e., it has no up-down symmetry) and stochastic. Surface 
inhomogeneities increase with time, signaling a non/mear evolution. However, since there 
are no significant overhangs, the growth mechanism should be essentially local; that is, 
the growth rate at a given point depends only on the surface pattern in the neighborhood 
of that point. (Kinetic roughening can also produce quite different surface patterns with 
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branched tree-like structures. Their growth is strongly nonlocal since the large trees shield 
the smaller ones from further deposition of material [Q.) 

Of course, generic scale invariance far from equilibrium is not limited to interfaces. 
Other important examples are hydrodynamic turbulence or slowly driven systems with 
so-called self-organized criticality: dynamical processes such as the stick-slip motion of 
an earthquake fault generate a power law distribution of "avalanches" with long-ranged 
correlations in space and time. A simple lattice model with a self-organized critical state 
is the so-called forrest fire model 0, ^. On a given lattice site, a tree grows with a small 
probability per unit time. With an even smaller probability, the tree is hit by a lightning 
which destroys it along with all the trees in the same contiguous forrest cluster These 
events are the avalanches. The dynamics leads to a self-similar stationary pattern of 
forrests and voids. 

A satisfactory theory of non-equilibrium scale invariance should classify the different 
universality classes and provide calculational methods to obtain the scaling exponents 
exactly or in a controlled approximation. To establish such a theory is obviously a complex 
task that will challenge statistical physicists probably over the next 25 years. The power 
law correlations on large scales of space and time should again be described by continuum 
field theories, although the proper continuum formulation is far from clear for many 
dynamical systems defined originally on a lattice. It is also an open issue which of the 
presently known field-theoretic concepts will continue to play an important role. For 
example, perturbative renormalization may fail to produce a fixed point describing the 
large- distance regime, as will be shown below for the example of driven surface growth. 
The seemingly easier task of describing the time-independent scaling in a stationary state 
is still involved since there is no simple Hamiltonian generating these correlations. 

This article collects a few results that may become part of an eventual field theory of 
nonequilibrium systems. We limit ourselves to models related to the Kardar-Parisi-Zhang 
(KPZ) equation || 



for a (i-dimensional height field h{r, t) driven by a force ri{r, t) random in space and time. 
This equation has come to fame as the "Ising model" of nonequilibrium physics. It is 
indeed the simplest equation capturing nevertheless the main determinants of the growth 
dynamics in Fig. 1(a): directedness, nonlinearity , stochasticity, and locality. The KPZ 
surface shown in Fig. 1(b) has been produced by a discretized version of the growth rule 
( |1.1|) and looks indeed qualitatively similar to these experimental data. 

The theoretical richness of the KPZ model is partly due to close relationships with 
other areas of statistical physics - notably hydrodynamic turbulence and disordered sys- 
tems - which are briefly reviewed below. Many more details can be found in refs. [0, ^. 
Due to the famous problem of quenched averages, disordered systems share some of the 
conceptual problems mentioned above. The observables are correlation functions aver- 
aged over the distribution of the random variables, which is not given by a Boltzmann 
weight. These correlations may be regarded as an abstract field theory but there is again 
no simple effective Hamiltonian. Such systems can have scale-invariant states at zero 
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temperature for which the very existence of a continuum hmit needs to be re-estabhshed. 

Despite considerable efforts, the KPZ equation has so far defied attempts at an exact 
solution or a systematic approximation in dimensions d > 1. The main reason is the 
failure of renormalized perturbation theory. Renormalization aspects of Eq. ( |1 . 1|) and 
its theoretical relatives will be discussed below. The main emphasis lies, however, on 
structures beyond perturbation theory. We analyze the internal consistency of the strong- 
coupling field theory expressed by the operator product expansion of its local fields. This 
approach turns out to be quite powerful. Using phenomenological constraints and the 
symmetries of the equation, it produces a quantization condition on the scaling indices 
from which their exact values in ci = 2 and d = 3 can be deduced. This quantization is 
somewhat reminiscent of what happens in conformal field theory and suggests the KPZ 
equation possesses an infinite-dimensional symmetry as well. 

It is not clear to what extent the results carry over to other nonequilibrium systems. 
Yet, the approach used here is fairly general and should be applicable in a wider context. 
It transpires that incorporating nonequilibrium phenomena into the framework of field 
theory will require yet another shift of focus to nonperturbative concepts and methods. 
This is likely to change our view of field theory as well. Negative scaling dimensions, 
dangerous variables, anomalies etc. are oddities today but may become an essential part 
of its future shape. 



1.1 Directed growth, Burgers equation, and polymers 

The time evolution of a KPZ surface depends only on the local configuration of the surface 
itself (and not, for example, on the bulk system beneath the surface). Hence, the r.h.s. of 
Eq. ( |1.1|) contains only terms that are invariant under translations h ^ h + const.: 

(a) The dissipation term V^/i is the divergence of a downhill current and acts to smoothen 
out the inhomogeneities of the height field. 

(b) The nonlinear term (Vh)"^ arises from expanding the tilt dependence of the local 
growth rate and acts to increase the inhomogeneities of the surface. A linear term b ■ Vh 
would be redundant since it could be absorbed into a tilt h ^ h + h-r. The higher powers 
{VhY, . . . turn out to be irrelevant in the presence of the quadratic term, as well as terms 
containing higher gradients such as (V^/i)^ or V^h. 

(c) The stochastic driving term 7]{r, t) describes the random adsorption of molecules onto 
the surface. It is taken to have a spatially uniform Gauss distribution with correlations 
only over microscopic distances. 



r/(r,t)=0, ri{r,t)ri{r',t') = a^6{r-r')6{t-t') . (1.2) 



A uniform average //(r, t) =rj would again be redundant since it could be absorbed into 
the transformation h{Y, t) h{r, t) — fjt. 

Eq. ( |1.1| ) is by no means the only model for a driven surface, and many experimental 
realizations of crystal growth are probably governed by related equations with additional 
symmetries and conservation laws [110] or with different correlations of the driving force 



see, for example, |[TT|). As the simplest nonlinear model, however, the KPZ equation 



remains a cornerstone for the theoretical understanding of stochastic growth. 
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The morphology of a rough surface is characterized by the asymptotic scahng of the 
spatio-temporal height correlations. In a stationary state, the mean square height differ- 
ence is expected to take the form 

mr^,h) - h{v^M)f) - |ri - r^r^g f , \ ; (1.3) 

\|ri - r2|^/ 

the higher moments {{h{Yi,ti) — /i(r2, ^2))^^) are of similar form. (In a system of finite 



size i?, Eq. ( |1.3|) is vahd for |ri — r2| <C -R and \ti — ^2! <^ ■) The scaling function Q 
parametrizes the crossover between the power laws {{h{ri,t) — h{r2,t)Y) ~ |ri — r2p^ and 
{(h{v,ti) — h{Y,t2)Y) ~ |ti — ^2^^^^^ of purely spatial and purely temporal correlations, 
respectively. These relations define the roughness exponent x > and the dynamic 
exponent z. In the marginal case x = 0; the surface may still be logarithmically rough. 



A surface governed by the linear dynamics (1J-) with A = has x = {2 — d)/2 and 
z = 2: it is rough for d = 1, marginally rough for d = 2 and smooth for d > 2. The 
phase diagram is well known also for A 7^ 0. In dimensions d < 2, any small nonlinearity 
(A/2)(V/i)^ is a relevant perturbation of the linear theory and induces a crossover to a 
different rough state called the strong-coupling regime. For d > 2, a small nonlinearity 
does not alter the smooth state of a linear surface. There is now a roughening transition 
to the strong coupling regime at finite values ±Ac IT^, |TB|, III . 



This phase diagram corresponds to the following renormalization group flow. For 
d < 2, the Gaussian fixed point (A = 0) is (infrared-)unstable, and there is a crossover 
to the stable strong-coupling fixed point. For d > 2, a. third fixed point exists, which 
represents the roughening transition. It is unstable and appears between the Gaussian 



fixed point and the strong-coupling fixed point, which are now both stable [|15], |T6|, |T^, |T8 

In one dimension, the critical indices of the strong-coupling regime take the exact 
values X = 1/2 and z = 3/2 [^9|, ^ Their values in higher dimensions as well as the 
properties of the roughening transition have been known only numerically p3| , 0, 

BB, n7\, pS[ and in various approximation schemes PPI, BUI, E^. In particular, it has 



been controversial whether there is a finite upper critical dimension c/> at and above which 
KPZ surfaces are only marginally rough (x = and z = 2). These issues are discussed in 
detail in Sections 3 and 4. 

Experiments on growing surfaces require delicacy since crossover and saturation effects 
can mask the asymptotic scaling. However, several experiments have produced scaling 
consistent with KPZ growth. The fire fronts in slow combustion of paper have x = 0.50 
and z = 1.53 in very good agreement with the KPZ values in d = 1. A recent study 
of kinetically roughened Fe/Au multilayers obtains x = 0.43 ± 0.05, which should be 
compared to the current numerical estimate x ~ 0.39 for d = 2 and to the presumably 
exact value ( [4.74| ). 

Eq. ( |1.1| ) is formally equivalent to Burgers' equation 

^t^r + (v ■ V) V = u V\ + V?7 (1.4) 

for the driven dynamics of the vortex-free velocity field v(r,t) = Vh{r,t) describing a 
randomly stirred fiuid (with A = —1) |]Ty|. In this formulation, Galilei invariance becomes 
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obvious: the substitutions 

h{T,t) h{r - ut,t) + u -r -^uH , v(r, t) ^ v(r - ut, t) + u (1.5) 

leave Eqs. ( |1 . 1|) and ( |1.4|) invariant. Due to this invariance, the roughness exponent and 
the dynamical exponent obey the scaling relation 



X + z = 2, (1.6) 

which guarantees that the total derivative dt = dt + u-V behaves consistently under scale 
transformations. 

It should be emphasized, however, that the velocity field of a stirred Burgers fiuid 
looks quite different from the gradient field of a KPZ surface since the driving force in a 
fiuid is correlated over macroscopic spatial distances R. This generates turbulence 



(somewhat different, of course, from Navier-Stokes turbulence). The velocity correlations 
show multiscaling. For example, the stationary moments ((f |j(ri) — f||(r2))'^) of the longi- 
tudinal velocity difference have a fc-dependent singular dependence on i? for |ri — <^ R. 
Multiscaling is not expected for driving forces with short-ranged correlations ( [1.2D . This 
point will become important below. 

Via the well-known Hopf-Cole transformation. 



Z(r, t) = exp 



2v 



;i-7) 



Eq. ( |1.1|) can be mapped onto the imaginary-time Schrodinger equation 



with 13 = l/2v 0. The solution can be represented as a path integral 

Z(r', t') = jvr 5{Y{t') - r') exp(-/3S) (1.9) 



with the action 



S= /*'dtf^(^)'-Ar^(r(t),t)| , (1.10) 



describing a string or directed polymer r(t) (i.e., the world line of a random walk) in 
the quenched random potential Xri{r,t) at temperature j3~^ = 2v. A configuration of the 
string is shown in Fig. 2. The stochastic driving term now appears as quenched disorder 
in a (1 + (i)-dimensional equilibrium system. This system is of conceptual importance as 
one of the simplest problems with quenched disorder. 
The transversal displacement of the string. 



((r(ti)-r(t2))2)~|ti-t2|'S (1.11) 

defines its roughness exponent C,. (Averages over the disorder are denoted by overbars, 
thermal averages by brackets (...).) The rough strong-coupling regime ( |1.3| ) of the growing 
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Fig. 2: A configuration r{t) of a string (or directed polymer) in a medium witli quenched point 
disorder. Due to tlie inliomogeneities of the medium, a typical path takes larger excursions to 
the left and to the right than an ordinary random walk. 



surface corresponds to a superdiffusive state of the string {( > 1/2) [g6], |38|. In this 
state, the universal part of its free energy in a system of longitudinal size L and transversal 
size R has the scaling form 

F(L, R) ~ L'^TiLR-^'^) . (1.12) 
In particular, the "Casimir" term 

7(i?) = lim dLF{L,R) ~ (1.13) 

L—*oo 

measures the free energy cost per unit of t of confining a long string to a tube of width 
R. The exponents ( and u: are related to the growth exponents by 

C = 1/z , uj = x/z. (1.14) 



The scaling relation ( |1.6|) now reads 

w = 2C-l. (1.15) 

Replica methods yield the exact exponents ( = 2/3 and uj = 1/3 for d = 1 but fail in 
higher dimensions. 

In the superdiffusive state, the free energy acquires an anomalous dimension —uj < 0. 
(At an ordinary critical point, the free energy is scale-invariant {u = 0), which implies a 
set of hyperscaling relations. Such relations are no longer valid for quenched averages.) 

The disorder-induced fluctuations ( p. . 1 1| ) persist in the limit f3~^ 0, that is, in the 



ensemble of minimum energy paths ro(t). In the weak-coupling (high-temperature) regime 
for d > 2, thermal fluctuations dominate {( = 1/2) and hyperscaling is preserved (cu = 0). 
The roughening transition between these two phases takes place at a finite temperature 
For d > dy, the Gaussian exponents govern the low-temperature phase as well, 
albeit with possible logarithmic corrections. 
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1.2 Overview of this article 



As emphasized already KPZ growth defines a field theory that is non-Lagrangian and non- 
perturbative. This article focuses on exact properties of its local correlation functions. 
We mention only briefiy the results of various approximation schemes. In particular, 
functional renormalization and references therein) and mode-coupling theory (|3l|, ^ 
and references therein) are important theoretical tools in a number of strong-coupling 
problems but their status in field theory is not yet fully understood. 

The first part of this article describes directed strings. These systems have a fascinating 
spectrum of physical applications and the language of directed strings proves to be an 
ideal framework to address some of the theoretical issues of directed growth. 

In Section 2, we discuss directed strings in thermal equilibrium without quenched dis- 
order. A single such string describes a free random walk and is thus generically rough. In- 
teractions of a single string with an external defect or mutual interactions between strings, 
however, can induce a localization transition destroying the long-ranged correlations of 
the rough state. (De-) localization phenomena are an essential feature of strings and sur- 
faces [|^; they will appear in several contexts throughout this article. We use renormal- 
ized perturbation theory to derive the phase diagram of directed strings with short-ranged 



interactions and the critical behavior at the (de-)localization transition |^T], ^ The 



results are, of course, well known and can be derived in various other ways. The approach 
used here stresses that the response of the system to perturbations probes the correlations 
in the unperturbed, rough state. It is based on the operator product expansion of the local 
interaction fields, a familiar concept in Lagrangian field theory (see, e.g., ref. |I|). 

In Section 3, this approach is extended to the field theory of directed strings in a 
random medium. In the replica formalism, a single such string is represented by a system 
of many strings without quenched disorder but with mutual interactions. The localized 
many-string state corresponds to the strong-coupling regime of the random system. Per- 
turbative renormalization of these interactions turns out to produce the exact scaling 
at the roughening transition for 2 < d < A but fails to describe the strong-coupling 



regime [0, 0. Insight can be gained by studying several strings in a random medium 
with additional direct interactions. These probe the disorder-induced correlations in the 
scale-invariant strong-coupling regime. The temperature becomes a dangerous irrele- 
vant variable at the strong-coupling fixed point; this is the field-theoretic fingerprint of 
quenched randomness. The direct interactions can be treated in renormalized perturba- 
tion theory about that fixed point |^ assuming the existence of an operator product 
expansion. We find again (de-)localization transitions, which are relevant for various 
applications. Their critical properties are given in terms of the single-string exponents. 
Comparing the effect of pair interactions in the strong-coupling phase and at the roughen- 
ing transition of a single string then shows that the single-string system - corresponding 
to the standard KPZ dynamics - has an upper critical dimension (i> < 4 

Section 4 returns to growing surfaces. The dynamical field theory of KPZ systems and 
its renormalization are discussed. Perturbative renormalization of the dynamic functional 
is compared to the string renormalization of Section 3 [|18|, and it is shown why pertur- 
bation theory fails for the strong-coupling regime in d > 1. However, the scaling in this 
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regime can be studied directly using the operator product expansion of the height field. 
We find that the KPZ equation can have only a discrete set of solutions distinguished 
by field-theoretic anomalies Comparing this set with numerical estimates of the 

exponents x ^i^nd z then gives their exact values in d = 2 and d = ?>. 



2 The field theory of directed strings 



Ensembles of interacting directed strings describe a surprising variety of statistical sys- 
tems in a unifying way. Examples are interfaces between different bulk phases in a 2D 
system steps on crystal surfaces |Q, fiux lines in a type-II superconductor or ID 
elastic media Directed strings are also related to mathematical algorithms detecting 
similarities between DNA sequences ||53|, |5^, ^ ^ ^ . 



At finite temperatures, a single string (or a collection of independent ones) would 
simply perform Gaussian fiuctuations. It is the interactions of the strings with each other 
and with external objects that generate the thermodynamic complexity of these systems. 
Attractive forces lead to wetting transitions of interfaces, bunching transitions of steps, 
and depinning transitions of fiux lines. All of these are transitions between a delocalized 
high-temperature state with unconstrained fiuctuations and a localized low-temperature 
state whose displacement fiuctuations are constrained to a finite width ^ ; for a review, see 
ref. EDI. In this and the next Section, we discuss a few such systems, emphasizing their 
common field-theoretic aspects. 

A single thermally fiuctuating string is given by the partition function 



Prexp(-/3S[r]) 



with the Gaussian action 



S[r] 



dt 



1 /dr^ 

2 I dt 



for the (i-component displacement field r(t). In a finite system (0 < t < L, < ri, 
i?), the universal part of the free energy has the scaling form 

F{L,R) =J^{L/(3R^) . 

This defines in particular the Casimir amplitude 

C{R) = (3^R^ lim ^^^(L, R) , 

L— >oo 



(2.16) 

(2.17) 
■ ,rd< 

(2.18) 
(2.19) 



measuring the scaled free energy cost per unit of t of confining a long string to a tube 
of width R. It depends only on the boundary conditions in transversal direction. For 
periodic boundary conditions, C = 0. 

The displacement field r{t) has the negative scaling dimension —(q = — 1/2. Its two- 
point function 

du —j^ (2.20) 
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Fig. 3: A thermally fluctuating directed string r{t) and a rigid linear defect at r = 0. Contact 
interactions between these objects are described by the local scaling field ^{t). 



requires an infrared regularization by appropriate boundary conditions. It is the difference 
correlation function 



((r(ti) - rit,)r) = -2(r(ti)r(t2)) + (At,)) + {r'it,)) ~ \h 



(2.21) 



that remains well-defined in the thermodynamic limit L, R oo and becomes scale- 
invariant. The exponent Co is called the thermal roughness exponent. 

Consider now a directed string interacting with a rigid linear defect at r = as shown 
in Fig. 3. If the interaction decays on the microscopic scale |r| ~ a, the system has the 
action 



(2.22) 



in the continuum limit a — 0. The local interaction of the string with the defect is 
proportional to the contact field $(t) = 6{r{t)) of canonical scaling dimension xq = d(o. 
The conjugate coupling constant go has the dimension yo = 1 — xq with t as the basic 
scale. 

Of course, this system can be treated exactly, for example by solving the imaginary- 
time Schrodinger equation 



/9- 



(2.23) 



for the wave function ( p..9| ); see |^ in the context of directed strings. Here we discuss a 



different way of solution ^ that can be generalized to problems with quenched 

disorder. For the purposes of this Section, it is convenient to set /5 = 1, which amounts 
to the substitution t — > (3~H in the action ( p.22| ). 



2.1 Correlation functions and the operator product expansion 

The perturbative analysis of the interaction in ( p.22| ) is based on the correlation functions 
($(ti) . . . $(tAr)) in the unperturbed state, which can be calculated exphcitly. We take 
each component of the displacement vector to be compactified on a circle of circumference 
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Fig. 4: Operator product expansion of contact fields for a thermally fluctuating string. The 
short-distance asymptotics of the pair of fields ^{t) and ^{t') is given by the single field 
times a singular prefactor. The dashed line indicates the string configurations generating the 
singularity \t — t'\~^°. 



R. The scale R also serves to generate the renormalization group flow defined below. 
With this regularization, longitudinal translation invariance emerges for "bulk" values 
-C ti, . . . ,t]\f -C L in the limit L — > oo independently of the boundary conditions at 
t = and t = L. 

The translation invariant one-point function 

{^{t)) = ($) = /?-^'«/fo (2.24) 

is simply the probability (density) of finding the fluctuating string r{t) at the origin r = 
for a given t. Similarly, the iV-point function ($(ti) . . . $(tAr)) is the joint probability of 
the configurations with N intersections of the origin at given values ti, . . . ,1^- 

These correlation functions develop singularities as some of the points approach each 
other. For example, the joint probability of intersecting the origin r = both at t and at 
t' equals the single-event probability ( p.24| ) times the probability of return to the origin 
after a "time" \t — t'\, which becomes singular as t' —>■ t: 

($(t)$(t')) = Co\t - ($(t)) + . . . (2.25) 

with Co = (27r)-^o. 

The structure of this singularity and the coefficient Co are local properties: they 
appear in any (connected) A^-point function as two of the arguments ti,tj approach each 
other, independently of the other points remaining at a finite distance and of the infrared 
regularization. This can be expressed by the field relation Hi 



$(t)$(t') = Co\t - m) + ... (2.26) 

illustrated in Fig. 4. It is called an operator product expansion (a familiar concept in 
field theory; see, e.g., ref. [||]). The dots denote less singular terms involving gradient 
fields. Such terms are generated, for example, if the r.h.s. of (|2.26| ) is expressed in terms 



of $(t') = $(t) + {t' — t)$'(t) + . . ., which leaves the leading singularity invariant. 
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2.2 Renormalization 



It is convenient to set up the perturbation theory for the Casimir amphtude ( |2.19| ). Since 
C is a dimensionless number, the contribution of the interaction 

ACK)=C((7o,i?)-C(0,i?) (2.27) 

depends only on the dimensionless coupling constant 

uo = goB^''^' . (2.28) 

The perturbation expansion 

AC(mo) = -R' E / . . . dt^mh) . . . ^t^)r . (2.29) 

contains integrals over connected correlation functions of the contact field in the Gaussian 
theory {qq = 0). Hence, the singularities of the operator product expansion (|2.26|) lead to 
poles in (^1291) . Inserting (|2]2|) and (|]26D into (|]2|), we obtain 

AC{uo) = ($) l^uo - %lj + 0{yX, 4) . (2.30) 

The same type of singularity (with different combinatoric factors) occurs in the expansion 
of correlation functions (^(ti) . . . $(tAr))(uo). For example. 



muo,R) = E^^/dti...dt^W)$(ti)...$(t^)r 

N=0 ■ 



= {<^)(^l-^uo^+0{y',uo,ul). (2.31) 

Perturbative renormalization consists in absorbing the singularities of the "bare" se- 
ries (|2.29|) and (|2.31|) into new variables up and $p defined order by order. (We use the 
subscript P to distinguish perturbatively defined couplings and fields from their nonper- 
turbatively renormalized counterparts. This distinction is not necessary in the present 
context but will become crucial below.) To leading (one- loop) order, the coupling constant 
renormalization can be read off directly from Eq. ( ^.30| ). Defining 

Up = ZpUQ (2.32) 

with 

Zpiup) = 1 - —Up + 0{y%p, uV) (2.33) 
the Casimir amplitude as function of up becomes finite to order Up, 

AC(up) = i?^o/fo(<l>) Up + O^y'^^ul.ul) . (2.34) 



12 



The coupling constant renormalization ( |2.33| ) implies a renormalization of the conjugate 
field, 

$p(t) = Zp^{t) (2.35) 



with 



Zpiup) = = 1 + — Mp + Oiy^up, Mp) . 
dup yo 



(2.36) 



This renders also the correlation functions (<l>p(ti) . . . ^p(tN)){up) finite. The scale de- 
pendence of Up is governed by the flow equation 



iip = CqROr Up 



youp 



1 — Mp(d/dMp) logZp 



Using ( |2.33| ), we obtain 



Up = yoUp - CoUp + 0{yoUp, Up 



(2.37) 



(2.38) 



Hence, the one-loop ^-factors and the resulting flow equation are entirely determined by 
the constants yo and Cq encoding local properties of the unperturbed theory. 

For this particular system, the one-loop equations turn out to be very powerful because 
the pole at order Uq is the only primitive singularity for y^ in the bare perturbation 
series. Hence, the theory is one-loop renormalizable [BDj , P^ : it can be described by 
the flow equation 

Up = yoUp — CoMp (2.39) 

terminating at second order. This property leads to exact results for local observables of 
the perturbed theory. In the Appendix, it is derived for the more general many-string 
system of Section 3.2. 

Of course, the form ( p.39| ) of the flow equation is not unique. It depends on the in- 
frared regularization of the bare perturbation series and on the renormalization conditions 
deflning the coupling constant up. Changing either amounts to flnite reparametrizations 



of Up. Linear reparametrizations change the coefficient of Wp in ( p. 39 ), while nonlinear 



reparametrizations introduce higher order terms (see the discussion in the Appendix). 
However, local observables of the perturbed theory are "gauge invariant"; i.e., indepen- 
dent of these choices. They can be computed exactly from (|2.39|) and the associated 
2-factors ( [A.13| ). The simplest such observable is the anomalous dimension 



X* = Xq — Up - — logZpiup) 
dup 



1 + 1/0 • 



(2.40) 



It governs the correlation functions of the contact field at the nontrivial fixed point Up of 
the fiow equation, for example, 

($(t)) ~ (2.41) 



and 



($(t)$(t')) ~ |t-tT"'($(t)) 



(2.42) 
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for |t — t'l ^ -R^. By simple scaling arguments, it follows from ( p.41| ) that the normalized 
string density 

P(|,'|) , <^W)-^')) (2.43) 



has the singularity 
with 



/dr(5(r(t)-rO) 
P(r) ~ r^R-'^-^ for r <^ R 
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X 

T 



d = Ayo = 2(2 - d) . 



(2.44) 



(2.45) 



Of course, the exponent 6 can be obtained in a simpler way. The Schrodinger equation 
has the singular ground state wave function Z ~ 7-2-^^ g^j^^ P Z^. 



2.3 Results and discussion 



For d < 2, the Gaussian fixed point = is unstable under the flow ( p. 391) and governs 



the (de-) localization transition. Close to the transition, the localization width has the 
singularity 

e ~ (-^7o)-^°/^° (go < 0) . (2.46) 

The fixed point u* is stable and determines the asymptotic scaling with repulsive contact 
interactions in the limit ^ 00 or (70 ^ 00. The string density P(r) then has a long- 
ranged depletion given by (|2.44| ) with 6 > 0. 



At the borderline dimension d = 2, where the two fixed points coalesce, the theory is 
asymptotically free with 

e ~ exp(-CoCoM) (^7o<0). (2.47) 
For 2 < d < 4, the transition is governed by the nontrivial fixed point with 

^--i9c-go)-'^"^'* {9o<9c) (2.48) 

and y* = 1 — X* = —yo- This fixed point now has a negative value of 6, describing a 
divergence of the string density P{r) for r — due to the attractive interaction. 

It is obvious that the same scaling occurs in a number of related systems. For a 
directed string (ri > 0,r2, . . . ,rd)(t) confined to half space, the term go6{ri) describes 
a short-ranged interaction with the boundary of the system. The scaling of the string 
close to the boundary is described by the above two fixed points for d = 1 with the 
transition point shifted to a value Qc < and the fluctuations parallel to the boundary 
decoupled. Similarly, a single directed string interacting with a linear defect is equivalent 
to the two-string system Z = J dridr2 exp(— /?S'[ri, r2]) with the action 

containing pair interactions \l'(t) = 6{ri(t) — r2(t)). The center-of-mass fluctuations are 
decoupled, and ( p.43| ) is the pair density with r(t) = ri(t) — r2(t). In d = 1, Gaussian 
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Fig. 5: (De-)localization in a system of many strings. In this example, the strings are steps on a 
crystal surface coupled by inverse-square and short-ranged forces. Typical step configurations at 
different temperatures: (a) Well above the critical temperature Tc, the steps are dominated by 
the no-crossing constraint and the long-ranged repulsion. Hence, they are well separated from 
each other with relatively small fluctuations, (b) In the critical regime near Tc, the probability 
of a step being close to one of its neighbors is substancially enhanced. This is accompanied by 
increased step fluctuations and a broader distribution of terrace widths, (c) Below the faceting 
temperature, the steps form local bundles. On average, the distance between two neighboring 
bundles is larger than the width of an individual bundle. The fluctuations of these "composite" 
steps are smaller than those of individual steps. 



strings with contact repulsions are known to behave asymptotically like free fermions |6T 



and ( |2.45| ) then gives the correct scaling P{r) ~ of the pair density imposed by the 
antisymmetry of the fermionic wave function. 

Letting aside these specifics of Gaussian strings, the qualitative features of the phase 
diagram are seen to rest on two properties of the local interaction field its scaling 

dimension is depends on ci in a continuous way and it obeys an operator product expansion 
with a self-coupling term ( p.26| ). These properties are found to be preserved for directed 
strings in a random medium despite the lack of a local action. 

It turns out that the renormalization discussed in this Section is also applicable to tem- 
perature-driven transitions in systems of many directed strings. An example of current 
experimental interest is vicinal surfaces, i.e., crystal surfaces miscut at a small angle with 
respect to one of the symmetry planes. A vicinal surface can be regarded as an ensemble 
of terraces and steps The steps are noncrossing (fermionic) directed strings. They 
are coupled by mutual forces that turn out to have a short-ranged attractive and a long- 
ranged repulsive part |62]. Typical step configurations are shown in Fig. 5. At high 



temperatures, the ensemble of steps is homogeneous. Below a faceting temperature (that 
depends on the step density), the steps are found to form local bundles. Hence, the 
surface splits up into domains of an increased and temperature-dependent step density 
alternating with step-free facets |Q . A critical regime containing the faceting transition 
(the analogue of the (de-)localization transition of two strings) separates the high- and 
low-temperature regimes. In the renormalization group, one still finds a pair of fixed 
points linked by an exact one-loop renormalization group for the strength of the contact 
interaction. These fixed points determine the faceting transition and the high-temperature 
regime, respectively The long-ranged forces influence the universal features (e.g., the 
scaling of the contact field) of both fixed points in a characteristic way. This determines, 
for example, the distribution of terrace widths and the density of steps in a bundle. The 
theoretical results compare favorably with recent experiments on Si surfaces 
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Fig. 6: Directed strings in a disordered medium with additional contact interactions, (a) A 
string and a rigid linear defect, (b) Two strings with mutual interactions, (c) A string and a 
wall. 

The thermodynamic complexity of this many-string system is due to an interplay of 
attractive interactions, Fermi statistics, and repulsive forces. In the next Section, we 
shall discuss the much simpler case of bosonic strings (i.e., strings allowed to intersect) 
with contact attractions only. In = 1, the latter will collapse to a bound state at any 
temperature. However, qualitatively different behavior emerges in the formal limit of 
vanishing number of strings. This limit turns out to describe a single string in a random 
medium. 



3 Directed strings in a random medium 



In this Section, interactions between strings play a double role. On the one hand, a single 
string in a quenched random medium can formally be represented as a system of p strings 
without disorder but with mutual interactions, in the limit p ^ [^. This well-known 
replica formalism turns out to be a convenient basis for the perturbative renormalization 
of the random system 



g. 



On the other hand, additional interactions are important in many applications of 
directed strings in random media. An example is the physics of superconductors [0, 
68| , |69| , [70| . A type-II superconductor in a magnetic field h above a critical strength 
hci contains magnetic flux lines at a density that depends on h. The lines are directed 
parallel to the magnetic field. Their thermally activated transversal fluctuations dissipate 
energy at the expense of the supercurrent. The superconductor may be doted with point 
impurities, linear or planar defects designed to suppress these fluctuations by localizing 
the flux lines. (Linear defects are generated, for example, by irradiating the material with 
heavy ions 0, 0].) 

As the external field approaches the critical value hd, the ensemble of flux lines be- 
comes dilute. It is then useful to study the approximation of a single flux line interacting 
with a single columnar defect in the presence of point impurities |]73|, |73|, ^ ^ ; see 
Fig. 6(a). The next step is to consider pair interactions in a dilute ensemble of lines as 
shown in Fig. 6(b) [^, [ 



80|, HI 



The interaction of strings with the boundaries of the system (Fig. 6(c)) can be dis- 
cussed on a similar theoretical footing. The particular case of one transversal dimension. 
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where the string becomes an interface of the system, has been of interest as a simple model 
for wetting in a random medium The interaction of a string with a linear defect is 
also relevant in the context of DNA pattern recognition |^ . Two related 

DNA sequences in different organisms have mutual correlations inherited from their com- 
mon ancestor in the evolution process. The algorithmic detection of these correlations 
turns out to correspond to a localized state of a string. 

It is not surprising that a random medium, by changing the displacement statistics of 
a single string, modifies also its direct interactions with other strings and with external 
objects. In the renormalization group for these transitions, quenched impurities enter in 
a characteristic way: the strong-coupling fixed point has a dangerous irrelevant coupling 
constant that alters the scaling properties of the direct interactions Consequently, 
the (de-)localization transitions differ from those in pure systems. In turn, the response 
to such interactions becomes a theoretical tool to study the correlations at the strong- 
coupling fixed point. This is used at the end of this Section to show that the single-string 
system has an upper critical dimension less or equal to four [E^. 



3.1 Replica perturbation theory 

A medium with quenched pointlike impurities exerts a local random potential ri{r, t) on a 
directed string. For a given configuration of the impurities, the partition function of the 
string is 

Z[ri] = J Vrexp{-pS[r,ri]) (3.50) 

with the action ( |1.1U| ), 

S[r,v] = +Vm,t)j , (3.51) 

where A has been set to 1. We take the local potential variables to have the Gaussian 
distribution given by (|1.2| ) and compute average quantities like the free energy 



F = -p-^ J Vr] exp J dtdr ^V^i^, r)) log . (3.52) 

At first sight, this system looks quite different from those of the previous Section. 
However, we can write logZ[?7] as a partition function of p independent strings labeled by 
the index a, 

logZy = lim- ( TT - 1 ) . (3.53) 

P^o P Va=l / 



In the replicated system, the integration over the rj variables can be performed |3^ . This 
gives 

F = lim - F„ (3.54) 
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with 

= log Jvr,... Vvp exp{-pS[v,, rj) (3.55) 



and the action 



S[ru...,r,] = /clt(|E^^)'-/5^'E'^«/^wj • (3-56) 

The coupling of the original string r(t) to the random medium now appears as a con- 
tact attraction ^afsit) = <^(i"a(t) — r^(t)) between the rephcas (i.e., phantom copies) 
Tiit), . . . , Tpit). The free energy of the rephcated system is related to the cumulant ex- 
pansion of the random free energy ||37[| , 

^p = E|!^^'^ (3-57) 

k=l 

The physical properties of this system strongly depend on p. For p > 1, the attractive 
interaction reduces the fluctuations of the strings. A single string {p = 1) undergoes no 
interaction. Randomness enhances the string fluctuations, and this is naturally associated 
with values p < 1. The existence of the random limit p — is not clear a priori. It can 
be established for d = 1, where the replica system is exactly solvable (see Section 3.2). In 
higher dimensions, the replica interaction can still be treated in perturbation theory. 

For arbitrary values of p, the path integral in Eq. ( p.55| ) can be rewritten in second 
quantization []TS|, 



Df^Dcj) exp 




/? / dtdr \(t)\dt- V 1 - 



(3.58) 



The contact attraction is described by the (normal-ordered) vertex 0^0^. Since this inter- 
action conserves the number of strings, the dependence of ( ^.5^ ) on p is contained entirely 
in the boundary conditions at early and late values of t. Hence, the boundary conditions 
are important for obtaining the replica limit p — > 0. 

The representation ( |3.58| ) is a convenient framework for perturbation theory [jl8[. We 



now mark the parameters /3(7^,(Tq, the free energy and all longitudinal lengths with the 
subscript to distinguish them from their renormalized counterparts introduced below. 
In the Appendix, the renormalization is carried out for the disorder-averaged Casimir 
amplitude 

C = PlR^J,{R) (3.59) 
defined by Eq. (|1.13| ). The disorder- induced part 

/SC{al (3^\ R) ^ C{al (3^\ R) - C(0, (3^\ R) (3.60) 

can be expanded in powers of the dimensionless coupling constant Uq = g^R^^^'^^'', where 

90 = -PX (3.61) 
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and yo = {2 — d)/2. Due to proliferation of replica indices, the perturbation series is more 
complicated than its analogue in the previous Section However, as shown in the 



Appendix, it is still one-loop renormalizable: the Casimir amplitude ( p. 591 ) is a regular 
function of the coupling constant up = ZpUQ defined by (|A.13|) , 

/3^AC(mp) = -^up + 0{youp, uj,) . (3.62) 

An immediate consequence of the one-loop renormalizability is that the strong-coupling 
regime is beyond the reach of the loopwise perturbation expansion ( |A.3|) since the flow 
equation ( p. 391 ) does not have a stable fixed point at negative values of |]18[. (The fixed 
point Mp = ?/o/Co > for (i < 2 is unphysical in this context since a repulsive interaction 
between replicas translates into a purely imaginary random potential.) We come back to 
this failure of perturbation theory in Section 4.2. 

The fixed point Up<0 for c? > 2 is to be identified with the roughening transition. 
At this fixed point, the Casimir amplitude (|3.59| ) takes a finite positive value without any 
explicit dependence on R, 

The scaling properties at the transition can hence be obtained exactly from the one- 



loop renormalization group [0, For example, the displacement fluctuations are still 
only diffusive, 

C = 1/2 , uj* = 0. (3.64) 
This follows by comparing (HD with the scaling C ~ R^^/'^ at a generic fixed point given 



by ( |1.13|) . Other exponents do depend on d. Conjugate to cTq is the local field 

$,(to) = J dr'r]{r', to)6{r{to) - r') , (3.65) 



which encodes the random potential evaluated along the string [0. Small variations of 
(Tg (i.e., of Up) are a relevant perturbation at the roughening transition. The dimension 

y = (3-66) 

is given by the eigenvalue of the beta function at the fixed point Up. The dimension of $^ 
is therefore x* = {4 — d)/2. In the action ( |3.56|) , this field becomes the replica pair field 
$Q,/3(t). The same dimension x* then follows from (|2.4CI| ) with the field renormalization 
( |A.13| ). One may also define the pair contact field \E'(t) = 6{ri{t) — r2(t)) of two real 



copies, i.e., two independent strings ri{t) and r2(t) in the same random potential. It can 
be shown that this field and its conjugate coupling constant also have dimensions x* and 
y*, respectively. 

The perturbative calculation can only be trusted for d < 4. The fact that x* would turn 
negative for d > 4 is clearly unphysical. Moreover, physical quantities become singular 



as d approaches 4; for example, C ~ ^/A — d M]. This shows that = 4 is a singular 



dimension for the roughening transition and ties in with the existence of an upper critical 
dimension (i> < 4 of the strong-coupling phase. 
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3.2 The strong coupling regime 

In the strong-coupling regime (i.e. for low temperatures or large disorder amplitudes), 
the string has the disorder-induced fluctuations 



((r(t) -r(t'))2) ~ (3.67) 

leading to anomalous scaling of the confinement free energy (|1.13|) , 

7(i?) ~ , (3.68) 

see refs. The exponents satisfy the scaling relation (|1.15| ). Superdiffusive scal- 

ing (tu = 2C — 1 > 0) is believed to persist up to an upper critical dimension dy (see 
Section 3.6). 

In d = 1, the exponents can be obtained exactly from the replica approach |3^. The 
system of p strings is always in a bound state for integer p > 1. The binding energy in a 
system of infinite width R, 

= lim dL (FpiPo, crl L) - Fp(/3o, 0, L)) , (3.69) 

can be computed by Bethe ansatz methods; one finds Ep ^ p + 0{p^). Analytically 
continued to p = and inserted in p.57| ), this gives F^ ~ L, hence u = 1/3, and x = 2/3 
by (OH)- 

In higher dimensions, the strings form a bound state only for > ctq^. If we assume 
Ep is still analytic in p and has the form Ep p + 0{p^°~^^) {ko = 2, 3, . . .), the same 
argument yields 

a; = ; (3.70) 

A;o + 1 

see the discussion in [|66| , |67[] . The exponents of the random system are indeed quantized. 



as will be discussed in the context of the Kardar-Parisi-Zhang equation in Section 4. The 
quantization condition ( [4.1| ) is consistent with ( |3.7U| ). 

In the continuum theory (|3.50|) , the large-distance scaling ( p.67|) and ( p.68| ) is reached 
in a crossover from diffusive behavior on smaller scales. The crossover has characteristic 
longitudinal and transversal scales to and f beyond which the disorder-induced fluctuations 
dominate over the thermal fluctuations. The dependence of these scales on the effective 
coupling ( |3.61| ), 



r 



'-go)-^/y' {d<2,go<0) 
(3^% = { exp(-Co/^?o) {d = 2,go< 0) (3.71) 
{9c -go)-'^''* {2<d<A,go<gc) 



with yo = 2/{2 — d) = —y*, can be obtained from the replica action ( |3.56| ), see ( |2.46D 
( p. 4^ ). The string displacement and the confinement free energy have the form 



((r(to) - r(t^))2) ~ 13^% n{to/to) (3.72) 
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and 

-2 D-2 



/o ~ (3^'R-'HRlr) , (3.73) 

with scaling functions that are finite in the short-distance hmits to and R <^ C,, 

respectively. In the opposite limit, comparison with ( p. 671 ) and ( |3.68| ) exhibits the singular 



dependence on the bare parameters Pq and ctq contained in the scaling functions, 

((r(to)-r(t^))2) ~ Po'%'nto - tX^ (3.74) 

and 

7o ~ P^'f-'R-' . (3.75) 
These singularities can be absorbed into the definition of the renormalized quantities 

t = {P/Po)to , 7 = {Po/P)Jo (3.76) 
written in terms of the renormalized temperature 



see the discussion in [^. The renormalized displacement function and confinement free 
energy remain finite in the continuum limit f — > (i.e. Pq^ — >■ or CTq — >■ oo), as follows 
by inserting ( |X75| ) and (^) into and (§TSp. 



The existence of a zero-temperature continuum limit is crucial if the ensemble of 
ground states generated by the quenched disorder is to have universal features. According 
to Eq. ( |3.77|) , the renormalized temperature is an irrelevant coupling constant of 
dimension —u. This is why the renormalized theory may be called a zero-temperature 
fixed point. It will be shown below that is a dangerous irrelevant variable which 



modifies the correlations of other fields in a characteristic way |45 . 



3.3 A string and a linear defect 

A single string coupled to a random medium and a linear defect has the partition function 
( p.5U| ) with the action 

S[r, v]=Jdt (^^ (^^) + vm,t) + (3.78) 

containing the contact field $(t) = 6{r{t)). The string configurations are determined 
by a competition between two kinds of interactions. Point defects roughen the string 
and make its displacement fiuctuations super diffusive. An attractive extended defect, 
on the other hand, suppresses these excursions and, if it is sufficiently strong, localizes 
the string to within a finite transversal distance ^. As in Section 2, the two regimes are 
separated by a second order phase transition where the localization length ^ diverges. In 
contrast to temperature-driven transitions, it involves the competition of two different 
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configuration energies rather than energy and entropy. Hence, the transition persists in 
the zero-temperature hmit. 

The properties of this zero-temperature phase transition have been controversial 



7^ , [75t [fq , Following the treatment of ref. we expand the defect contribution 

/\C{g,R) = C{g,R)~C{Q,R) (3.79) 

to the zero-temperature Casimir amplitude 

C = r(^-^)/<J{R) (3.80) 

about the point g = given by the strong-coupling continuum theory. This leads to a 
perturbation series formally analogous to ( |2.29| ), 

AC{g, R) = Y: / dt2 . . . dt^(<l>(ti) . . . $(t^))^ . (3.81) 

Of course, the perturbation series cannot be evaluated explicitly since the multipoint 
correlation functions of the contact field at the zero-temperature fixed point are not 
known exactly. However, one can still write down the one-point function 

($(t)y = J^= R-'^'^ (3.82) 

(with X = dQ and establish the short- distance structure of the higher connected correla- 
tion functions. 



Consider first the displacement function ( |1 . 1 1| ) . It has the low-temperature expansion 

((r(t) -r(t'))2) = \t- tf ^ + [3-^ \t - + . . . , (3.83) 

assuming analyticity of the crossover scaling form ( p.72 ) in the scaling variable 
Hence at zero temperature, ( |1.11| ) equals its thermally disconnected part (r(t) — r(t')y. 
The connected part can be shown to equal that of the Gaussian theory ||38[| , 



{{r{t)-r{t')yyr.r'\t-t'\; (3.84) 



it appears as the leading correction to scaling in (p.83|). 



In analogy to ( |2.25 ), the two-point function ($(t)$(t')) is assumed to factorize for 
\t — t'\ <^ R^/'^ into the one-point function ($(t)) times the i?-independent return prob- 
ability to the origin, which is proportional to the inverse r.m.s. displacement given by 
(S): 

($(t)<i>(t)) ~ |t - t'r"(i - (3~^d'^\t - t'r" + . . .)($(t)) . (3.85) 

Again the leading singularity is due to sample-to-sample fluctuations of the minimal en- 
ergy paths, while the correction term is due to thermal fluctuations around these paths. At 
zero temperature, the fleld $(t) can be replaced by its thermal expectation value ($(t)); 
hence ($(t)$(t')) equals its thermally disconnected part ($(t)) (<I>(t')) and the connected 
part {^{t)^{t')Y vanishes, just as the connected displacement function (|3.84|) does. The 
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,t) 



Fig. 7: Operator product expansion of contact fields for a string in a disordered medium. The 
short-distance asymptotics of the pair of fields ^{t) and ^{t') is given by the single field <I>(t) 
times a singular prefactor. The dashed lines indicate the string configurations generating the 
singularity \t — t'\~^~^. 



subleading term in ( p.85| ) is the sum of ^{t')y and a temperature-dependent correc- 
tion to ($(t)) (^(t')). An analogous argument applies to the singularities in any correlation 

0. This leads to the operator product expansion 



function (. . . $(t)<l>(t') . . .) as |t - t'\ 



with a constant C > 0. Defining the contact field $(r',t) = S{r(t) 
generalized to the spatio-temporal operator product expansion 



<l>(r,t)<l>(r',t') r^Cp-^\t-f 



/ I —X — LJ 



'u\t - t'l 



<l>(r,t) 



producing a spatial singularity of the form 

$(r,t)$(r',t) ~/3->-r'|-^-"/^<l>(r,t) . 



(3.86) 
( |3.86| ) can be 

(3.87) 
(3.88) 



The operator product expansion encodes the statistics of rare fluctuations around the path 
of minimal energy |7^ as shown in Fig. 7. Notice that the leading singular term is 
proportional to the irrelevant variable I3~^ and hence governed by a correction-to-scaling 
exponent. This is why the temperature is called a dangerous irrelevant variable. 

As before, the operator product expansion ( p.86| ) determines the leading ultraviolet 
singularities of the integrals in ( 3.81 ), which appear as poles ixiy = 1 —x — uj. The defect 
contribution to the Casimir amplitude can again be written in terms of a dimensionless 
coupling constant: 



with 



u 



u u 

K y 



(3.89) 



(3.90) 

Indeed, Eq. ( p.89|) remains finite in the limit (3^^ for fixed R and g. The pole 
can be absorbed into the "minimally subtracted" coupling constant um = Zmu (which 
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should not be confused with the couphng ( |2.32| ) defined in perturbation theory about the 
Gaussian fixed point). With 

C 

Zm{um) = 1 um + 0{u\i) , (3.91) 

y 

the fiow equation iiM = CR^rUm reads 

iiM = yuM - Cull + 0{ul^) . (3.92) 

We conclude that an attractive linear defect in a random system is less effective in 
localizing a directed string than in a pure system, in agreement with some previous ap- 
proximate renormalization group studies |7^, |7^, ^ . A weak defect defect is a relevant 
perturbation of the zero-temperature fixed point um = ^ only in dimensions d < 1. It 
localizes the string with 

e ~ {-g)-^''' (^7 < 0) . (3.93) 
At the borderline dimension = 1, the theory is again asymptotically free, with 

e ~ eM-CCla) (^7 < 0) . (3.94) 

For d > 1, a finite defect strength \g\ > \gc\ is necessary to localize the string: 

e~(^7c-^7)-^/^* ig<gc), (3.95) 

where = hm„^,^^„*^ um/ — um) = —y + 0{y^) is the eigenvalue of the fiow equation 
at the nontrivial fixed point. The disorder- averaged string density ( |2.43D has the short- 
distance singularity P(r) ~ with 9 = 2y/( + 0{y'^) < 0. 

A weak repulsive defect is an irrelevant perturbation of the strong coupling fixed point 
for d > 1. In particular, it does not generate a "fermionic" zero of the string density P(r) 
in the limit R —>■ oo. This prediction of renormalized perturbation theory is consistent 
with the scaling at a strong repulsive defect or barrier, which can be obtained exactly 



in d = 1 ||8^. An impenetrable barrier {g oo) is equivalent to a system boundary, 
leading to a long-ranged depletion P(r) ~ r of the string density (see Section 3.4). A 
hardly penetrable barrier has rare crossings that happen whenever the difference in the 
random energies on the left and right sides of the barrier in a given longitudinal interval At 
exceeds the barrier penalty. Since the string configurations one side consist of essentially 
uncorrected pieces of length ~ i?^/^, this difference scales as 

AF(At, R) = (At)i/3_^d(Atp-3/2) ~ (At)i/2i?-i/4 (3.95) 

for At ^ R^^"^. Hence the path remains on one side of the barrier for a typical longitudinal 
distance At given by (At)^/^i?~^/^ ~ g. This determines the expectation value of the 
contact field, 

M ~ ^ ~ g-'R-'^' . (3.97) 

We conclude that the penetrability g~^ conjugate to $ is a relevant perturbation with 
eigenvalue j/ = 1/3 at the barrier fixed point g~^ = 0. It induces a crossover to the 
delocalization fixed point g = 0: the barrier becomes irrelevant on large scales. 
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3.4 A string and a wall 

As mentioned above, a half-space Gaussian string (ri > 0,r2, . . . ,Td)(t) in contact inter- 
action with a system boundary at ri = is equivalent to the string ri(t) in full space 
coupled to a defect at ri = 0. In a random background, this is no longer the case. The 
disorder potential couples the string coordinates ri, . . . ,rd and the boundary has a non- 
local influence on the string by cutting off all disorder configurations in the half space 
ri < 0. Alternatively, the half-space system can be understood as an unrestricted system 
with a defect plane at ri = and "mirror" constraints r]{ri, . . . ,rci,t) = ri{—ri, . . . ,rd,t) 



on the random potential |83| 



We restrict ourselves here to the case d = 1, where a half-space string with the action 
( p. 7^ ) can be treated exactly by Bethe ansatz methods or by mapping on a lattice 



gas H^. One finds a (de-)localization transition with the localization length singularity 

^--{9c-9)-' {9<9c). (3.98) 



At the transition, the disorder-averaged string density (|2.43| ) has the singularity ( |]40 
p. 317) 

P(r) ~ r-^/^i?"^/^ for f < r < P. (3.99) 

Hence, the boundary contact field $b(t) = S{r(t)) has the expectation value 

{<l>b{t)) ~ R-'^"/^ (3.100) 

with Xb = 1/3. By arguments as in Section 3.3, we then obtain the operator product 
expansion 

Mt)Mt') = Cbf3-'\t - t'l-^'^-^^bit) + ... (3.101) 

with Cb > 0. Hence, the wall changes the statistics of rare fluctuations. This operator 
product expansion leads again to a beta function of the form (|3.92| ) with i/b = I — Xb — uj. 

In d = 1, the interaction with the wall is a truly relevant perturbation with eigenvalue 
?/6 = 1/3 at the (de-)localization point. The leads to a singularity ( |3.93|) of the localization 
length, which agrees with (|3.98| ) and with the result of obtained from a variational 
scaling argument for the bound-state free energy. For g > Qc, there is a crossover to 
the fixed point u\.j = Ub/Cb + Oiyl) with exponents x* = 1 — oj + Ub + Oiyl) and 6 = 
x*/( — d > describing a long-ranged depletion of the string density P(r). This is again 



in agreement with exact results and numerical transfer matrix studies but the 



one- loop calculation underestimates the true value 6 = 1. 

3.5 Strings with mutual interactions 

The effective action 

1 / dr- ^ 



S[r„r,,r^] = Jdtl-j:[^) +J:vir^{t),t)+gm) (3-102) 



1=1 \ / 1=1 
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with the contact field \l/(t) = S{ri(t) — r2(t)) describes a pair of directed strings that five 



in the same random medium and are coupled by short-ranged mutual forces |177| , [Tq , |f9| , 



5D|, In the case of fiux lines, for example, this interaction is repulsive. Due to the 
random potential, the center-of-mass displacement of the strings does not decouple from 
their relative displacement. The two- string system is therefore not equivalent to a single 
string and a linear defect. 

The effects of pair interactions turn out to be much stronger than in a pure system. 
Qualitatively, this is easy to understand. At zero temperature, two noninteracting strings 
in the same random potential share a common minimal path ro(t). At small but finite 
temperatures, it turns out that the strings still follow a "tube" of width f around ro(t) 
with finite probability. Hence, their overlap probability (^^(t)) remains finite in the limit 
i? — s> oo, in contrast to that of noninteracting thermal strings, (^'(t)) ~ R~^'^". This 
explains the strong sensitivity of the system to repulsive forces. However, the strings 
make large individual excursions from the tube. The disorder-averaged pair density ( |2.43| ) 
(with r = ri — has the form [^1], ^ 



P(r) ~ for r > f and ^ oo (3.103) 

dictated by the operator product expansion ( |3.88| ). A strong repulsion {g oo) forces 
one of the strings onto the lowest excited path ri(t) that has no overlap with ro(t) (with 
fiuctuations of the form ( |3.103D around ri(t) at finite temperatures). The paths ri(t) and 



ro(t) have an average distance of order R, and the pair density should have a long-ranged 
depletion P(r) ~ with 6 > for g oo, just as that of free fermions in d = 1. At 
the strong-coupling fixed point of noninteracting strings, the pair field \E'(t) is therefore a 
relevant perturbation inducing the crossover to the "fermionic" behavior for g ^ oo. 
This argument can be made quantitative |H6|. We start from an expansion 



AC{g, R) = f; fdh... dt^ i^mih) . . . ^{t^)r (3.104) 

N=l ■ •' 

of the Casimir amplitude ( p. 801) about the strong-coupling fixed point of noninteracting 
strings [g = Q). Since the overlap probability of the two strings 

(3.105) 

remains finite as i? ^ oo, the local field \E'(t) has dimension x = 0. The connected 
correlation functions of \E'(t) vanish at zero temperature, as do those of the contact field 
$(t) of Section 3.3. One obtains an operator product expansion of the form ( p.86| ), 

^(t)^(t') = C2p-^\t - t'|""^(t) + . . . . (3.106) 

with C2 > 0. Its leading term is again a correction to scaling proportional to the dangerous 
irrelevant variable f3~^. Renormalization of the perturbation series (|3.104|) then leads to 
the fiow equation ( p.92| ) with y = \ — oj. 

It follows that an attractive pair interaction always localizes the strings. The pair 
density P(r, of the localized state has the form 

P{r,i) r-. (3-\-<^~^/'^V{r/i) forr>f; (3.107) 
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the scaling function V has a finite hmit at short distances and decays exponentially on 
scales r^^. The localization length ^ has the singularity 

e ~ {-9)-^'' = < 0) . (3.108) 

The transversal scale defined by the mth moment of the pair density, 

/ /• \ l/m 

^ / dr r- P(r, (m = 1, 2, . . .) , (3.109) 



scales as 

^ /5-i/"^(_^)-{C-W'")/?/ (^ < 0) . (3.110) 

Recall that at an ordinary fixed point, all the scales have the same exponent as 
the correlation length ^. The dangerous irrelevant variable (3~^ breaks this universality 
and induces the multiscaling ( 3.110] ). A similar phenomenon occurs for thermal directed 



strings in dimensions d > 4 ||58| . 

Repulsive forces lead to an asymptotic scaling P(r) ~ with 9 = x*/( — d given in 
terms of the dimension x* = 2(1 — uj) + 0((1 — ujY) of the pair field at the nontrivial 
fixed point. A numerical transfer matrix study of a pair of strongly repulsive strings 



{g oo) confirms the long-ranged depletion of the pair density [85|. The exponent 9^2 
is underestimated by the one-loop calculation. 

The above considerations are valid only as long as d is below the upper critical di- 
mension dy of a single string. This can be seen as follows. It is possible to show that 
the ground state path ro(t) is unique (up to microscopic degeneracies of the order of the 



lattice spacing) in almost all realizations of the disorder [^8|, This uniqueness is also 
manifest in the pair density (|3.103|) : for any fixed ro > 0, the probability of finding the 
strings at a relative distance r > tq remains finite for i? — oo, and that limit value tends 



to zero for /3 ^ — ^ |^ , 



dr P(r) ~/3-Vo"'^/'^ . (3.111) 

/r>ro 

However, for d ^ d^ (i.e., a; ^ 0), the pair distribution ( p.l03| ) shows a singular broad- 
ening: for P — >■ cx) and fixed the probability (|3.111 ) approaches one. Consequently, 



the overlap probability ( p.l05| ) goes to zero, {'^{t)) ~ u. This suggests that the statistics 
of string configurations becomes more complicated for d > d^. The strings no longer 
cluster in the vicinity of the minimal path as expressed by ( |3.111|) , but exploit multiple 
near-minimal paths at any finite temperature. Their overlap (^^(t)) is expected to vanish 
for R —!■ oo. 



3.6 Upper critical dimension of a single string 

The theory of strings with pair interactions discussed in Section 3.5 has an important 
implication for the single-string system |^ : the upper critical dimension of a single string 
is less or equal to four. As d^ is approached from below, the exponents ( and u tend to 
the Gaussian values ( = 1/2 and u = continuously. Hence, the upper critical dimension 
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could serve as the starting point for a controlled expansion. The name "upper critical 
dimension" is, however, quite misleading since (i> does not mark the borderline to simple 
mean-field behavior as in the standard theory of critical phenomena. On the contrary, the 
state of the system in high dimensions may even be more complicated, having presumed 
glassy characteristics 0. 

The existence of an upper critical dimension has been controversial. Numerical work 
seems to indicate that a strong coupling phase with nontrivial exponents 2; < 2, x > 
persists in dimensions d = A and higher ||2^, However, the results for d > 3 are 

not very reliable since the available system sizes are limited and corrections to scaling are 
not taken into account 

Various theoretical arguments favor the existence of a finite upper critical dimension dy 
at or slightly below four. Most of these approaches rest on approximation schemes (such as 
functional renormalization p9|, ^ or mode-coupling theory |^T|) whose status is not very 



well understood. The same is true for a recent approximate real-space renormalization [ P7[ ] 
challenging these results. 

The argument given here |^ is of a different nature; it is not tied to any of these 
approximation schemes. An important ingredient is ( 3.110 ), a set of exact relations in the 
strong-coupling regime. These relations describe the bound state (|3.107| ) of attractively 
coupled strings in terms of the only independent single-string exponent u = 2( — 1. We 
focus on the temperature dependence of the scales For fixed g, is monotoni- 

cally increasing with temperature according to ( |3.110| ). This is not surprising since it is 
temperature-driven fiuctuations of the strings around the path ro(t) that generate the pair 
distribution ( |3.107| ). At the roughening transition (i.e., for = the singularity of 
^^"^^ changes, 

^ ^ ^ i-gy^^^y* with y* = id- 2)/2 , (3.112) 

as discussed in Section 3.1. With the natural assumption that remains a mono- 
tonic function of for all P'^ < [3^^ and fixed g, one then obtains the inequalities 
y/ ((^ — uj/m) > y*/C*. These imply an upper bound on the free energy exponent: 

4- d , , 

(3.113) 



d 



The result dy < 4 then follows immediately. 

It is tempting to speculate about the nature of the strong-coupling regime in high 
dimensions. Below dy, the pair distribution of noninteracting strings at fixed temperature 
has the finite limit ( p.l03| ) for R — >■ 00, and this limit distribution collapses to S{ri — 
for — > 0. For d > d^, the pair distribution is expected to depend on both and 
R in an essential way. Its asymptotic behavior will then depend on the order in which 
the zero-temperature limit and the thermodynamic limit R 00 are taken. Similar 
properties are familiar from glassy systems. 



3.7 Discussion 

Interacting strings in a random background turn out to have a rich scaling behavior. In one 
transversal dimension alone, there are no less than six universality classes characterized by 
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of the exponent 9 of the disorder-averaged string density and the renormahzation group 
eigenvalue y of the contact couphng. These universahty classes describe 

• the (de-)localization transition of a string at a linear defect (6' = 0,?/ = 0), 

• the (de-) localization transition of a string at an attractive wall {9 = —l/2,y = 1/3), 

• the (de-)localization transition of a pair of strings with contact attraction {9 = 
-3/2,y = 2/3), 

• the scaling of a string at a barrier [9 = l,y = 1/3), 

• the scaling of a string at a repulsive wall {9 = l,y = —2/3), 

• the scaling of a pair of strings with contact repulsion {9 = 2,y = —4/3). 

Clearly, the list could be continued by including higher dimensions, long-ranged in- 
teractions or disorder correlations etc. Recall from Section 2 that without quenched 
disorder, the six cases are described by just two universality classes, namely Gaussian 
strings {9 = 0,y = 1/2) and free fermions {9 = 2,y = —1/2). 

This scenario is consistent with an operator product expansion (|3.86|) of the contact 
field, which is the basis for perturbation theory about the strong-coupling fixed point. 
The existence of an operator product expansion is of conceptual importance since the 
correlations in the strong-coupling regime are generated by global minimization of the 
free energy and not by a local action. The operator product expansion explicitly contains 
the temperature as dangerous irrelevant variable. This variable, which is proportional 
to the surface tension u of the associated KPZ surface, will prove dangerous in the growth 
problem as well: it generates the dynamical anomaly discussed in the next Section. 



4 Directed growth 

All of the methods and results on directed polymers in a random medium discussed in 
the previous Section have their correspondences in KPZ surface growth via the Hopf-Cole 



transformation (|1.7| ). In particular, the KPZ equation has an upper critical dimension 
less or equal to four (see Section 3.6). 

To discuss dynamical renormahzation, we rewrite Eq. (|1.1|) as a path integral for 
the height field h{r,t) and the response field h{r,t). It proves necessary to distinguish 
carefully between renormalized fields and couplings defined in a nonperturbative way, 
and perturbatively renormalized quantities defined, for example, by minimal subtraction. 
The perturbative renormahzation of the dynamic path integral is seen to be equivalent 
to the replica perturbation theory of Section 3: it captures the roughening transition but 
does not produce a fixed point corresponding to the strong-coupling regime |T8|. This 



failure of perturbation theory is explained by the fact that the nonperturbatively renor- 
malized coupling constant u and the corresponding height field h (defined by conditions 
on correlation functions at a renormahzation point) have a singular dependence on their 
perturbatively subtracted counterparts up and hp. 

The strong-coupling regime thus calls for nonperturbative methods. In the second 
part of this Section, we analyze constraints on the effective growth dynamics imposed by 
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the consistency of correlation functions. The consistency relations take again the form 
of an operator product expansion. The response field h is found to obey an operator 
product expansion similar to those of the previous Section [EH]. The resulting structure 



of the response functions is tied to KPZ growth with additional deterministic driving 
forces. In particular, a surface with a local variation —{g/X)S{r) in the deposition rate 
corresponds to a directed string interacting with a linear defect (see Section 3.3). It turns 
out that (for A > 0) a sufficiently enhanced deposition {g/X < Qc/X) induces an increased 
growth rate even in the thermodynamic limit R —>■ oo, while for g/X > gc/X the growth 



rate is asymptotically independent of g [|^, |73[. These two regimes are separated by a 



nonequilibrium phase transition, the analogue of the string (de-) localization transition. 

Finally, we turn to the correlations of the height field ||4^. These are assumed to obey 
an operator product expansion as well. Consistently with the available numerical data, 
we further assume these correlations do not show multiscaling, unlike the correlations 
in a turbulent fluid. This property constrains severely the possible form of the operator 
product expansion. A discrete set of values of the roughness exponent emerges: 

X = -^ withA;o = l,2,.... (4.1) 
ko + 2 

Only for these values, the KPZ equation admits solutions with a dynamical anomaly in 
the strong-coupling regime. It is argued on phenomenological grounds that the scaling of 
growing surfaces should be governed by such solutions with an odd value of ko. Hence the 
exact values of the growth exponents for two- and three-dimensional surfaces are derived. 



4.1 Field theory of the Kardar-Parisi-Zhang model 



The path integral formulation of the stochastic evolution (|1.1| ), ( |1.2D is obtained in a 
straightforward way. We rewrite the Gaussian distribution of the driving force ?7(r, t). 



X'?7exp 



drdt 



2a2 



'Drj'Dh exp 



drdt 



h + hrj 



(4.2) 



introducing the purely imaginary "ghost" fleld /i(r, t). Eliminating rj by using the equation 
of motion, we obtain the (Ito-dicretized) path integral 



VKDh exp 



(4.3) 



the generating functional of the dynamical correlations. Ghost fleld insertions produce 
the response functions 



/ N _ N+N \ N 

nMri,t,) n Mr,>t,)) = n 



5 



N+N 

n K^v^i 

j=N+l 



where p(r, t) is an additional source fleld in the equation of motion, 

A 



dth = z/V^/i + ^(V/i)^ + V + P 



(4.4) 



(4.5) 
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There is a third way to express the KPZ dynamics. Writing the equal-time height 
correlations in the form 

(/i(ri, t) . . . /i(r„, t)) = (/i(ri) . . . /i(r„))t = Jvh /i(ri) . . . /i(r„) Pt , (4.6) 

the time dependence has been shifted from the field variables h{r,t) to the configuration 
weight Pt[h]. The latter obeys the functional Fokker-Planck equation 

where 

J(r) = uV^h{r) + ^{Vh)\r) (4.8) 

is the deterministic part of the current. 

In the linear theory, the fields h{r, t) and h{r, t) have canonical dimensions — xo = 
{d — 2)/2 and Xo + d = {d + 2)/2, respectively, and the dynamical exponent is zq = 2 (the 
basic scale is now that of r). 

A crucial role in the sequel will be played by the infrared divergencies of the height 
correlations Their physical meaning can be seen already in the linear theory, which 
has the response propagator 

wMr.t.)) - (J{;;:;;i)./. exp {-^^) (4.9) 

{0{t) denoting the step function) and the height-height correlation function 

(/i(ri, h)h{r2, h)) = a^J drdt {h{r, t)/i(ri, ti)) (/i(r, t)h{r2, h)) . (4.10) 

However, the last expression is divergent in infinite space for d < 2 (i.e., xo > 0). There- 
fore, it depends strongly on the initial or boundary conditions. For stationary growth in 
a system of size R, the dominant part in the limit i? — > oo is 

(/i(ri, ti)/i(r2, t2))R ~ {h^)R ~ i?2x ^ (4.11) 

where (/i^) = {h'^{ri,ti)) = {h'^{r2,t2)) with periodic boundary conditions. The amplitude 
( [4. lip measures the global roughness of the surface, i.e., the size of its mountains and 
valleys. The regularized height correlation, however, can be written in terms of height 
differences measuring the local roughness, 

(Mri,ti)/i(r2,t2))i? - {h^)R = -^((Mri,ti) - Mr2,t2))')R • (4.12) 
The finite limit for i? — > oo, 

((Mri,ti)-Mr2,t2))') = |ri-r2p^°^off^^^l , (4.13) 



(ri - r2)2^ 
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shows (affine) scale invariance. In the interacting theory, the infrared regularization is 
more comphcated since also the higher connected height correlations develop i?-dependent 
singularities. But the local statistics of the surface, measured by height differences, is still 
expected to remain finite for R ^ oo, i.e., to decouple from the divergent amplitudes. 

The interacting theory has an ultraviolet singularity as well. Taking the expectation 
value of the equation of motion, we obtain the growth rate, which has the form 

{dth)R = ^{{Vhy-)n ~ a'^-' + 0{R'^'') (4.14) 

for stationary growth. The dependence on the short-distance cutoff a can be removed by 
the subtraction 

(V/i)2 ^ (V/i)2 - ((V/i)2)oo , (4.15) 

which affects the one-point function {h{r,t)) but leaves the higher connected correlation 
functions invariant. This normal-ordering of the interaction field will always be implied in 
the sequel. The nature of the cutoff dependent term in ( [4.14| ) is clear from the mapping 



onto directed strings. Via Eq. (|1.7| ), the subtracted growth rate is directly related to the 
Casimir term ( |3.8(JD , 

-mh)B =7iR) , (4.16) 
while the a-dependent term contributes to the nonuniversal part of the free energy. 



4.2 Dynamical perturbation theory 



Perturbation theory for the KPZ dynamics can be set up in different ways: 

(a) The standard formalism is a diagrammatic expansion of the path integral ( ^.3| ) with 
interaction vertex h{WhY about the Gaussian limit A = 0. The propagators of the 
Gaussian theory are the response function {h{ri,ti)h{r2,t2)) and the correlation func- 
tion {h{ri,ti)h{r2,t2)); see Eqs. ( ^4.9| ) and ( [4.10| ). Renormalization is usually based 
on the expansion of Fourier transformed two-point functions {h{—'k,—uj)h{k,uj)) and 
{h{—k, —u!)h(k,u!)). This has been carried out to first [0, |T^ and second order [0, |T^; 
see also the criticism of ref. [0. In the Appendix, we discuss instead the expansion for the 
real-space response function (/i(ri, ti)/i(r2, ^2)) and the universal part of the growth rate 
{dth)R in order to exhibit the analogy to the replica perturbation theory of Section 3.1. 
Beyond leading order, however, calculations become cumbersome in this formalism. 

(b) By the Hopf-Cole transformation ( |1.7|) and the corresponding transformation of the 
response field 

/i = ^log0, h=^(p(p, (4.17) 



the path integral 



A 2u 
3|) can be brought to the form [^] 



V(j)V(f) exp 




~2i7 



12/2 



(4.18) 



which equals the replica action ( |3.58| ) introduced previously in ref. (with 13 = l/2u 
and A = 1). We emphasize, however, that ( [4. 18] ) is not a faithful representation of the 
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growth dynamics unless it is defined with boundary conditions enforcing the hmit p ^ 0; 
see the discussion in the Appendix. The theory is now Gaussian for finite A but cr^ = 0. 
The propagator (0(ri, ti)0(r2, t2)) is again given by (|4.9D; the interaction vertex 0^0^ 
arises from the noise term h"^. The resulting diagrammatic expansion is identical to the 
replica perturbation theory of Section 3 and can hence be renormalized exactly to all 
orders, as shown in the Appendix. 

In both cases, the expansion parameter is again the dimensionless coupling constant 
uo = goR^yo with 

^ U ''''^ 

where the fields and couplings of the bare theory are labeled by the subscript 0. The poles 
in Ho can be absorbed into loopwise subtracted fields hp, hp and the coupling constant 
Up = ZpUQ given by ( [A.13| ). There is no renormalization of to = ^p? which would be 
required in the strong coupling regime according to ( p.7(: ). 

The flow equation is again ( |2.39| ) with the fixed point u*p = yo/Co representing the 
roughening transition for 2 < d < 4. At the transition, the Casimir amplitude of the 
growth rate takes a finite value regular in i/q, 

R' J^(9,,hp)n = + 0{yl) , (4.20) 

which is directly related to the Casimir amplitude ( p.63|) of a string in a random medium. 
The scaling exponents at the transition, 

= , = 2 , (4.21) 

are independent of d, in agreement with the finite-order calculations of refs. |15, [16|, 
These exponents have been predicted previously by a scaling argument p3|. 



4.3 Renormalization beyond perturbation theory 

Of course, renormalization is not tied to a loopwise expansion and can be defined in 
a nonperturbative way [jl8|. The correlation functions show crossover scaling from the 
linear to the strong-coupling regime {d < 2) and from the critical to the strong-coupling 
regime (c? > 2). The crossover has the characteristic scales f and i given by ( |3.71|) with 
= 2u. 

For d < 2, the response and correlation functions have the scaling form 

/ N N+N \ 



Y[ho{rj,toj) n hoi^j^toj] 
^i=i j=N+i 



NN 



— Tfc toj ~ tok 



R ' R'o 

Eq. (|4.22| ) can be written as the bare Callan-Symanzik equation 



Uo 



R 



d_ 
dR 



d 



AO 3 



N 



d d 

N+N 



du, 



Nxo + N{xo + d) 







\i=i j=N+i 



(4.22) 



(4.23) 
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with 

Mo = RdR uo = 2yoUo . (4.24) 

For the infrared-regularized correlators, the explicit dependence on R vanishes in the 
thermodynamic limit. 

In the strong-coupling regime, these correlation functions develop anomalous scaling 
and hence a singular dependence on the bare coupling constant Uq. Renormalization 
consists in absorbing these singularities into new variables 



h 
t 
u 



Zhho 
Ztto , 
Zuq = 



h 



(4.25) 



The 2^-factors can all be written as functions of u. 

Under this change of variables, Eq. (^4.23|) transforms into the renormalized Callan- 
Symanzik equation 




Ao 



Nxiu) + Nixiu) + d) 



X 



\j=i 



j=N+l 



with 



u = ROr u 
z{u) 



2you 



1 — u{d/du) \ogZ 

zq - u—\ogZt , 
au 



du 



(4.26) 

(4.27) 

(4.28) 
(4.29) 



Notice that only two of the 2^-factors are independent. The ghost field renormalization is 
tied to that of h since the response function {h{ri, ti)h{r2, ^2)) always has dimension d by 
its definition. Furthermore, Galilei invariance implies the relation Z = ZJ^'^Z^'^ obtained 
by inserting the exponent relation x + ^ = 2 into ( [4.27|) , ([4.28|) , and ([4.29|) . A different 
but equivalent Callan-Symanzik equation is derived in ref. |[T^. 

The renormalized variables ([4.25|) can be defined in a nonperturbative way by two 
independent normalization conditions. These are imposed, e.g., on the stationary two- 
point functions in an infinite system. 



R^^°{{h{0,t)-h{R,t)f){u) 
R'^{h{r,t)h{r,t + R^)){u) 



i?2*'((/^o(0,to)-/io(^,to))')(0), (4.30) 
R''{hoiT,to)hoir,to + R'))iO) . (4.31) 



R is now an arbitrary normalization scale. An alternative to (|4.3(]| ) is the normalization 
condition 

R'-'''{dth)Riu) = -^u (4.32) 
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on the universal finite-size correction to the stationary growth velocity in a system of 
size R. Comparing ( |4.30| ) and ( [4.31j ) with the asymptotic scaling of the bare functions 



obtained from ( |4.22| ), 

((/io(ri,to)-/io(r2,to))') ~ f-2(™)|ri - r2p^ (4.33) 

(^o(r,toi)/io(r,to2)) ~ if'''\t,2-hir"\ (4.34) 
we infer the asymptotic behavior of the ^-factors: 

Zh ~ (f/i?)^-^" , Zt ~ {f/Ry-^^ . (4.35) 

The renormalized height correlation and response functions have a finite continuum limit 
f — >■ 0. Of course, this does not imply that all observables are finite in this limit since com- 
posite fields may generate additional singularities. However, as discussed in Section 4.5, 
the structure of composite fields of the KPZ theory is expected to be quite simple. All 
negative-dimensional fields (for example, the monomials h^ir)) have scaling dimensions 
given by the linear spectrum ( |4.51| ). Consequently, the transformations ( [4. 251 ) do renor- 
malize the correlations of these fields. 

The renormalized height field h{r,t) obeys again an equation of motion of the form 



( |1 . 1|) but the coefficients become singular in the limit r ^ 0. By substituting the renor- 



malized variables ( [4.25 ) into ( [L.l| ), ( [4. 3] ) or ( [4.7| ), we obtain 



Z^^uo ~ z/* X {f/RY , 
a' = ZlZ-'alc^a*'x{r/R)''-'+'^^ , (4.36) 
u = uqZ-'^Z^'^ ~ u* 

with finite limit values z/*, o"*^, and g* . Galilei invariance is refiected in the asymptotic 
scale invariance of the dimensionless coupling, m ~ m*. This is the nonperturbative 
strong-coupling fixed point. The other coefficients become irrelevant a.s f / R 0. In 
particular, we recover the scaling (|3.77|) of the renormalized temperature = 2z/. It is 
easy to verify that the strong-coupling asymptotics ( [4.35| ) and ( |4.36| ) remains unchanged 



for d > 2 although the scaling functions ( [4.22| ) now contain singularities marking the 
roughening transition. 

We emphasize again that the existence of this finite continuum limit depends only 
on the existence of a well-defined scaling regime at large distances. We have not as- 
sumed perturbative renormalizability, i.e., that the 2^-factors can be computed as power 
series in uq. It is instructive, however, to compare the renormalized theory with the 
perturbatively finite theory of Section 4.2. Recall that in a usual e-expansion, imposing 
normalization conditions like ( [4.30 ), ( 4.31[ ) is equivalent to requiring finiteness order by 



order in perturbation theory. The respective coupling constants, u and up, are related by 
a diffeomorphism that remains regular in the limit e —>■ and is defined on a domain of 
interaction space that contains the fixed points up = u = and u*{u*p). This equivalence 
is lost in the crossover to the strong-coupling fixed point. To see this, we integrate the 
flow equation ( p.3£ ) for up in two dimensions with the initial condition up{Rq) = Ui. The 



solution 

1 - ui\og{R/Ro) 
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diverges at a finite value R = Ri = Roexp{l/ui). It follows immediately that up{R) is 
not a good parametrization of the crossover. The pole of Up{R) is only a "coordinate 
singularity" of perturbative subtraction schemes; any renormalized coupling u{R) 
is expected to remain regular ai R = Ri. Hence the function up{u) has a pole at the 
value u{Ri) between the fixed points m = and u* . Accordingly, the loopwise subtracted 
correlation functions at the normalization point become singular as well: perturbative 
subtraction fails to maintain a smooth crossover of the large-distance regime. 



4.4 Response functions in the strong coupling regime 

The correspondence between directed strings in a random medium and growing surfaces 
extends to strings with additional interactions. We restrict ourselves here to the case 
of a string and a linear defect discussed in Section 3.3 to infer the properties of the 
renormalized ghost field correlations at the strong-coupling fixed point. 
The action (|3.78| ) translates into the growth equation 

^^h = uy^h + -{yhf + ri-^5{v) , (4.38) 
2 A 

see refs. [p8| , The extra term describes a spatial inhomogeneity in the average rate of 
mass deposition onto the surface. The Casimir amplitude (|3.79|) of the free energy is that 
of the growth rate, 

KC{g, R) = -\R'^+^ mh)R{g) - (9^/1)^(0)) . (4.39) 

For gf/A < (c? < 1) and g/X < gd^ {d > 1), translational invariance is strongly bro- 
ken. The string bound state corresponds to a surface state with a macroscopic stationary 
mountain 

H{r) = {hiO,t))-{h{r,t)) (4.40) 

generating an enhanced growth rate even in the limit R oo. 

The inhomogeneity in Eq. ( [4. 381 ) contributes a term — {g/X) J dt h{0, t) to the dynamic 
action in ( [4.3|) . Hence, the perturbation theory ( p.81| ) for the Casimir amplitude (|3.79|) 
can be reproduced in the dynamical formalism p5| , 

AC{g, R) = \ hm {exp[~{g/ \)Jdt h{0, t)] h{r, L)) . (4.41) 

_L— >oo 

The two series are related term by term via the mapping (|4.17|) [^ , 

(/?A)™-i m,)---^{tN)Y = lim (/i(0, ti) . . . /i(0, t^) hiv, L)) . (4.42) 

L^oo 

Eq. ( p.87| ) translates into an operator product expansion for the ghost fields, 

hir, t)h{r', t') = xc\t - n (^^rr^) ^) ' ^^-^2) 

which produces an equal-time singularity h{r,t)h{r' ,t) ~ |r — r'|~*^'^+'^) h(v,t). It is quite 
plausible intuitively that the response to a pair of nearby sources reduces to the response 
to a single source times a divergent factor. This form of the operator product expansion 
follows also from the mode-coupling approach of ref. 
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4.5 Height correlations in the strong coupUng regime 

In the sequel, we discuss the renormahzed and connected equal-time height correlations 

(/i(ri, t)... /i(r„, t)) = (/i(ri) . . . /i(r„))i (4.44) 

in the strong-coupling regime, following ref. An infinite surface growing from a 

fiat initial state h{r, 0) = develops height correlations depending on the differences 
Tij = r, — Tj and the correlation length ~ t^/^ increasing with time. (In a finite system 
of size R, the correlation length will eventually saturate at a value ^ ~ i?. We restrict 
ourselves here to unsaturated growth, i.e., to values R ^ ^j.) In the scahng regime 

r<|r,,|<et (i,j = l,...,n) , (4.45) 

the height correlations will generically become singular as some of the points approach 
each other. For d < d^, these singularities should be power laws. They are assumed to 
follow from an operator product expansion 

h{r,) . . . h{r,) = Y: \r^,\-'^'^+^- f • • • . ^1 (^(^i) , (4-46) 
o \\^12\ \ri2\J 

an identity expressing any ra-point function as a sum of [n — k+ l)-point functions in the 
limit \rij\ <^ |rj/| -C = 1, . . . ,k and I = k + 1, . . . ,n). The sum on the r.h.s. runs 

over all local scaling fields 0{r). Each term contains a dimensionless scaling function 
(a simple number for k = 2) and a power of |ri2| given by the scaling dimensions xo and 
Xh = —X (such that the overall dimension equals that of the l.h.s.). The field Ok with the 
smallest dimension, Xk, determines in particular the asymptotic behavior of the fc-point 
functions as — > oo, 

{hir^)...hivk))tr-{Ok)tr-^r' ■ (4.47) 

The amplitudes {Ok)t = {Ok{j))t diverge with ^t, i-e., Xk < 0. They measure the global 
roughness, which increases as the surface develops higher mountains and deeper valleys. 
Local properties of the surface should, however, behave quite differently. For example, 
the gradient correlation functions are assumed to have a finite limit 

lim (V/i(ri) . . . Wh{vn))t = (V/i(ri) . . . V/i(r„)) . (4.48) 

By writing h(ri) — h{r'^ = J^^ ds- V/i(s), the same property follows for the height difference 
correlation functions 

((Mn) - h{r[)) . . . (Mr„) - hir'J)), . (4.49) 

This implies a feature familiar from simulations: one cannot recognize the value of C,t from 
snapshots of the surface in a region much smaller than 

The stationarity condition ([4.48|) restricts severely the form of the height correla- 



tions ( [4.44| ). This can be seen as follows. The operator product expansion ( [4.46| ) induces 



an expansion for the gradient field v = V/i of the form 

v(rO . . . v(r,) = E |ri2|-^-+^'- C'f f • • • , 0(ri) (4.50) 

n V ri2 I"i2 / 
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Fig. 8: Operator product expansion for a /c-tuple of KPZ height fields. The hnes indicate 
the contraction of the fields, i.e., the coupling to the identity 1. Partial contractions would 
not contribute to connected correlation functions. Subleading singular and regular terms are 
omitted. 



with new scaling functions and the dimension = — x + 1. (Both sides of (^4.501) 
are tensors of rank k whose indices are suppressed.) The fields O on the r.h.s. govern the 
time-dependent amplitudes (v(ri) . . .\{Tk))t ~ {0)t ~ ^t~^^ in analogy to ( |4.471 ). Hence, 
the stationarity condition ( ^.48| ) allows in ( [4.50| ) only fields O with a non-negative scaling 
dimension xo, such as 1 (the identity field), (V/i)^(r), etc. This in turn restricts the 
possible terms in ( [4.461 ): 



(a) Singular terms involving fields 0{r) with xo > 0. In particular, we call the coefficient 
function of 1 the contraction of the fields h{ri), . . . , h{rk). 

(b) Regular terms, where the coefficient \ri2\~^^''~^^° is a tensor of rank in the 
differences rn {i = 2, . . . , k). Such terms do not violate (|4.48|) since they have a vanishing 
coefficient in ( [4.50| ) for N < k. They can readily be associated with composite fields 
of dimensions 

Xk,N = -kx + N . (4.51) 

The leading (A^ = 0) term involves the (normal-ordered) field Ok{r) = /i'^(r) and governs 
the asymptotic singularity ( 4.47| ) with Xk = x^fl = —kx- The higher terms correspond to 



fields with k factors h(r) and powers of V. 

The leading terms of ( [4.46| ) are shown in Fig. 8. It is useful to introduce the (normal- 



ordered) vertex fields 

Z,(r)=exp[g/i(r)], (4.52) 

which are the generating functions of the fields h^{r). Eq. ( [4. 461 ) then implies the operator 
product expansion 

Z,,iv,)Z,,iv,) = exp (^Ciy,w'}j Z,,+,,(ri) + OiC^f') , (4.53) 

where C^i = C^j(0, . . . , 0, ri2/|ri2|, . . . , ri2/|ri2|) with the first k arguments equal to 
and Wi = gi|ri2|'^ [^. This is a generalization of Wick's theorem. In an expansion of 
the exponential, each term Cl 1W1W2 represents a contraction of k + l fields h. Subleading 
singular terms (with positive-dimensional fields O) and regular terms (with fields con- 
taining height gradients) are omitted. Eq. ([4.53|) determines the asymptotic behavior of 
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the vertex n-point functions, 



[Z^,{r,)...Z^M)t-^expi^fY.<li] ■ (4-54) 



i=l 



If J2i Qi = 0, they have a finite hmit (Zg^(ri) . . . Z_q^..._g,^_^(r„)). Since these are precisely 
the vertex correlators that generate the height difference correlation functions ( 4.49D and 
since the vertex operator product expansion is analytic in the g^, this leads back to the 
stationarity condition ( ^.481 ). 

The operator product expansions ([4.46|) and ([4.53|) with the linear dimensions ([4.51|) 
are at the heart of the field theory for KPZ systems. They describe the following structure 
of height correlations: 

(a) The single-point amplitudes 

{h% ~ (4.55) 
are the moments of a probability distribution, 

{h')t = J dhh'^Pi^tih) . (4.56) 

The time dependence of the distribution takes the form 

A,(/i) = er^Pi(/^er'^) (4.57) 

expressing global scaling of unsaturated growth. Recall that by normal-ordering Eq. ( |1 . 1| ) , 
we have eliminated the nonuniversal part of the average height, {h)t ~ t. 

(b) Since ( |4.46| ) does not have singular terms with negative- dimensional operators, the 
local correlation functions ( [4.48|) and (|4.49|) decouple from the global amplitudes (|4.55| ). In 
particular, powers of the height difference hi2 = h{ri) — h{r2) reach stationary expectation 
values for — >■ oo determined by the contraction term in ([4.46|) , 

{hl^) ~ Irial'^^ . (4.58) 
These can be written as moments of a stationary probability distribution, 

{hl2) = /d/ii2/it2i^2(/ii2,|ri2|), (4.59) 

which has the scaling form 

P2{hu,r) = r-^V2{hi2r-'') . (4.60) 

Eq. ( [4.58| ) can be verified exactly for (i = 1 and seems to be consistent with the presently 
available numerical data also in higher dimensions. 

It is instructive to compare this structure with models of turbulence. Burgers' equation 
( |1.4|) with force correlations 

r/(r, t)r]{r', t') = erl6{t - A K ~ ^ J (4.61) 
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over large distances ro develops multiscaling: for example, the longitudinal velocity dif- 
ference moments 

(h(ri) -t;||(r2)]'=) ~ Irur'^--"'' r^'^ (4.62) 
have a fc-dependent singular dependence on |ri2| and ro for 

ro/7^< |ri2| <ro , (4.63) 

where TZ = r^'^e^l'^ jv denotes the Reynolds number |3^. This so-called inertial scaling 



regime should be compared to the KPZ strong-coupling regime f -C |ri2| -C ^t- According 
to ( |4.62|) , local and global scaling properties are no longer decoupled. Similar multiscahng 
is present in Navier-Stokes turbulence. Kolmogorov's famous argument predicts the exact 
scaling dimension of the velocity field, = —1/3, from dimensional analysis This 



determines the scaling of the third moment in ( [4.62| ) since = 0. The higher exponents 
X4,X5,... < cannot be obtained from dimensional analysis. Assuming the existence 
of an operator product expansion (|4.5CI|) , the term ( [4.62|) is generated by the lowest- 
dimensional field Ok with a singular coefficient. A similar operator product expansion is 



discussed in ref. Multiscaling thus implies the existence of a (presumably infinite) 

number of composite fields with anomalous negative dimensions. For the velocity vertex 
fields exp[gf (r)] of Burgers turbulence in one dimension, Polyakov has conjectured an 
operator product expansion similar to ( [4.53D P5| . The distinguishing feature of KPZ 
surfaces is the absence of multiscaling expressed by (f4.48|) . 



4.6 Dynamical anomaly and quantized scaling 

The operator product expansion ( |4.53| ) and the dispersion relation (4.51) have to be 
consistent with the underlying dynamics ( [4.7| ). However, as explained in ref. |^ for 
Burgers turbulence, the equation of motion for the renormalized correlation functions is 
quite subtle due to anomalies dictated by the operator product expansion. To exhibit the 
anomalies for the height correlations we introduce the smeared vertex fields 

Z^{r) = exp (^q J dr' 5„(r - r')h{r')^ (4.64) 

(where Sa{r) is a normalized function with support in the sphere |r| < a) and the abbre- 
viations = Z^^{r,), Z^ = Zg^(r,). 

Using ( [4.6| ) and ( [4.7[ ), it is straightforward to derive the equation of motion 

n 

{Zl . . . = Q^{Z1 ...JZt... ZD, , (4.65) 



1=1 

where 



JZt = [g.a^5,(0) + J(r,)]Zf . (4.66) 
The singularity structure of the current is determined by ( [4.46] ) and ( [4.36| ): 



oo 



JZt = g*Z, + a"^-^ [Y,^ Cka%^j Z, + 0{a^) (4.67) 
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for a, r — s> with a/f kept constant. The field Zq{Y) = (V/i)^Zg(r) denotes the finite part 
of the operator product (Vh)'^{r)Z^{r) for a — > 0, and Zi = Zg.(rj). The finite dissipation 
term (V^/i)Zg.(rj) becomes irrelevant in this limit since u ~ a^. The singular part of 
( [4.67| ) is a power series in Qi with asymptotically constant coefficients Ci = a*'^a'^6a{0) + 
^*ci,i + g*C2,i and = iy*Ci^k + g*C2,k for k = 2, 3, . . .. The terms of order a^^^''^^"^ 
originate from contractions V'^h{ri)h{r[) . . . /i(r^) ~ 1 and (VhY{Ti)h{r[) . . . /i(r'^) ~ 1. 
(Their respective coefficients ci^k and C2,k are complicated integrals involving the scaling 
functions in (|4.46| ), the regularizing functions Sa{ri — r'), and the ratio a/r.) Of course, 



divergent terms have to cancel so that Eq. ([4.65|) has a finite continuum limit 



dt {Zi... Zn)t = Y.(li(Zi...JZi... Zn)t (4.68) 

with 



JZ, = lim JZ^ . (4.69) 
For generic values of this implies jTZj = g*Zi. However, if x satisfies the condition 



( [4.1|) for some integer k^, the dissipation current contributes an anomaly: 

JZ, = g*Zi + u*ci,koqi''Zi . (4.70) 



Eqs. (|4.68|) and ( [4.70|) govern in particular the stationary state of the surface. For 



d = 1, the stationary height distribution P[h] is known, 

P ~ exp 1^-— J dr (Vh)^^ . (4.71) 

It equals that of the linear theory, thus restoring the up-down symmetry h{r) — {h)t — *■ 
—h{r) + {h)t broken by the nonlinear term in ( |1 . 1|) . The exponent x = 1/2 satisfies ([4.1|) 
with ko = 2 but the up-down symmetry forces the anomaly to vanish (ci,2 = 0). In higher 
dimensions, this symmetry is expected to remain broken in the stationary regime. The 
surface has rounded hilltops and steep valleys, just like the upper side of a cumulus cloud 
(argued in ref. [|100|| to be a KPZ surface). Hence, the local slope of the surface is correlated 
with the relative height, resulting in nonzero odd moments ((V/i)^(ri) [/i(ri) — h{r2)]''). 
However, this is consistent with Eqs. ( [4.6(^ ) and ( f4.70| ) only for odd values of fco, where 

(Z,(ri)^-,(r2)) - (Z_,(ri)Z,(r2)) = -{u* / g*) c,,,,q'^' {Z,ir^)Z^,ir2)) (4.72) 

and hence for odd values oi k > ko 

((V/^)^(rO[Mri) - h{r2)]') = -{v^ I g^) ci,.o([Mri) - Mr2)]'=-'«) • (4.73) 



The directedness of the stationary growth pattern thus requires a nonzero anomaly 



ko 



with an odd integer ko. The roughness exponent is then determined by Eq. (11). The 
values fco = 3 for (i = 2 and ko = 5 for d = 3 give the exponents ( [4. 741 ) and ( [4.75| ) 
quoted below. These appear to be consistent with the numerical results x ~ 0-39 and 



X ~ 0.31 p5|, |2^, respectively, and with the experimental value x = 0-43 ± 0.05 for 

d = 2 ra. 
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4.7 Discussion 



The scaling of strongly driven surfaces has been determined by requiring consistency of 
the effective large-distance field theory subject to a few phenomenological constraints. 
These are the existence of a local operator product expansion (4.46) and of a stationary 



state ([4.48|) that is directed (i.e., has no up-down symmetry). The stationarity condition 



has an important consequence: the decoupling of local and global scaling properties. The 
latter can be expressed by the time-dependent height probability distribution at a single 
point, the former by the stationary distribution of the height difference between nearby 
points. 

The Galilei invariance of the dynamic equation conspires with these constraints to 
allow only discrete values of the roughness exponent in two and three dimensions: x = 
2/(A;o -l- 2) with an odd integer k^. Comparison with numerical estimates then gives the 
exact values 

X = 2/^, 2 = 8/5 forrf = 2 (4.74) 

and 

X = 2/7, z = 12/7 ior d = 3. (4.75) 

The underlying solutions of the KPZ equation are distinguished by a dynamical 
anomaly in the strong-coupling regime: the dissipation term contributes a finite part 
to the effective equation of motion ( [4.68|) despite being formally irrelevant. The anomaly 
manifests itself in identities like ( [4.73| ) between stationary correlation functions, which 



can be tested numerically. 

There are at least three directions of future research where these concepts and methods 
can be of use. As mentioned in the introduction, the field of stochastic growth is rather 
diverse, and the theoretical cousins of the KPZ equation could be analyzed in a similar 
way, aiming at a better understanding of these non-equilibrium universality classes. In 
contrast to the KPZ equation, some of these systems do show multiscaling |101| , which 
makes their field theoretic description certainly more complex. 

Further applications to the theory of disordered systems are equally important. The 
present theory should be extendable from a directed string to higher- dimensional mani- 
folds in a random medium. 

Finally, this theory highlights both the theoretical similarities and differences of surface 
growth and fluid turbulence. While it is quite surprising that the scaling exponents of an 
interacting dynamical field theory can be predicted exactly in two and three dimensions, 
there is one other such case in a related system: the exact Kolmogorov scaling of the third 
velocity difference moment in Navier-Stokes turbulence. The higher velocity difference 
moments show multiscaling (see Section 4.5). For the simpler case of Burgers turbulence. 



dynamical anomalies are intrinsically connected to the multiscaling exponents [|102|] . This 
link is expected to be important in a wider context of turbulence. 

Perhaps some of these fascinating scaling phenomena far from equilibrium are not as 
inaccessible as they have appeared so far. 



42 



Appendix 



This Appendix contains some details about the renormahzation of the rephcated string 
system (Section 3.1) and the dynamic functional (Section 4.2). 
We start from the replica partition function 



VcpDcj) exp 



Lo 



dtn 



R 



dr Ufa 



+ 



2 J,2 



(A.l) 



with the effective coupling constant ( p.61| ), obtained from ( |3.58| ) by the change of variables 



^0 ~^ ^to- 
imaginary 

To evaluate 
define 



The field (j) is real- valued with the constraint > 0; the field is purely 
The normal-ordered interaction go(f)^(f)'^ conserves the number of strings. 
|A.1|) , we take periodic boundary conditions for the components of r and 



Zp = {p,r\Z\p) 



(A.2) 



with the initial state \p) = {l/p\)[J dr0(r)]^|O) at to = and the final state {p,r\ = (O|0^(r) 
at to = Lq. Correlation functions {■ ■ ■)p are defined in a similar way. With an appropriate 
normalization of the propagator, these boundary conditions normalize Zp{go, Lq = 0,R) = 
Zp{gQ = 0,Lq,R) = 1 for any value of p. 

The Casimir amplitude is defined in analogy to ( 2.19| ) and ( 2.27 ). Its expansion in 
terms of the dimensionless coupling uq = goR^y°, 



ACr, 



-go 



N=l 



J dto2 • • • dt0iv($2(0)$2(t02) • • • $2(t07v))p , (A.3) 



is a sum involving connected pair field correlations in the p-string sector of the unperturbed 
theory {uq = 0). The integrals in Eq. ( A.3 ) are infrared-regularized by the system width R; 
their ultraviolet singularities are determined by the short-distance structure of the pair 
field correlations and have to be absorbed into the coupling constant renormahzation. 
Hence consider the asymptotic scaling of the A^-point function ($2(^01)) • • • ^2(toN))p as 
the points toi,...,toiy approach each other. More precisely, we define to and Tj^ by 



to Tjk and let 



to/R^ 







(A.4) 



with Tjk and the "center of mass" t'o = N Hj=i toj remaining fixed. The asymptotic 
scaling is given by the p-independent operator product expansion 



N+l 



$2(t0l) . . . ^toN) = E to^""-'^^'^"/' [C^in, . . . , T^.,) $„(t[,) + . . .] 



(A.5) 



m=2 



where are scaling functions of the N — 2 linearly independent distance ratios tj^, 
the dots denote subleading terms down by positive integer powers of to/R"^, and we have 
defined the normal-ordered m-string contact fields 



$m(to) 



dr I 



r, tn 



(r,tn 



(A.6) 
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(a) (b) 



X3< XXX )CX 



Fig. 9: Singular diagrams contributing to the finite-size free energy ACp(uo) in tlie expansion 



( A.3 ). Tlie lines denote unperturbed single-line propagators (i^(ri, ioi)0(r2) ^02))) the dots pair 
contact vertices $2(^0)- (a) One-loop diagram containing a primitive pole in i/q. (b) A^-loop 
diagram containing a pole in yo of order — 1. 



For m = N = 2, ( |A.5| ) reduces to ( |2.26D ; i.e., C| = Cq. The operator product expansion 



dictates the leading singularities of the integrals in 

N-2 



[ dtot^-^ n ($2(^01) . . . ^2{toN))p = 
■' 1=1 

N+l 

Jn C+^^^'^-^+^^dto {^^{Q), + ... (A.7) 



m=2 

with 



. N-2 



Jn= f U dnC'^{n,...,TN_,). (A.8) 

1=1 

The diagrammatic representation of the perturbation series (|A.3| ) is discussed in detail 
in ref . 1 33 ] . At any integer value of p, the series is somewhat simplified since 



($m)p = for m > p. (A. 9) 

For p noninteger, all values of m contribute, e.g., to ([A.7|) . Hence, the series becomes 



complicated in the random limit p ^ 0. Its pole structure at yo = 0, however, remains 
simple. Consider the term in ( [A.5| ) with N = m = 2, corresponding to the diagram 



of Fig. 9(a). The loop has the value -R^^''c(?/o)/z/o with c{yo) = Co(?/o) + 0{yo). The 
pole originates from the universal short-distance singularity in ( [A.7| ), while the finite part 
depends also on the infrared regularization. Hence, we obtain to one-loop order 




with i?'^($2)p = p(p- l)/4. The pole can be absorbed into the definition up = ZpUQ with 
the .E-factor (|2.33|) , which leads to the beta function ( p.38|) and the field renormalization 
= Zp^ with (g). 



It is not difficult to discuss higher orders. The ultraviolet singularities of the A^-th 
order integral ([A.7| ) are contained in the coefficients or arise from the integration over 
to- In the first case, they are due to a proper suhdiagram and hence already absorbed 
into the renormalized coupling constant at lower order. Only the divergences from the 
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integration over to, with J"^ denoting the regular part of ( |A.8| ), may contribute to the 
primitive singularity at order A^. Inspection of ( |A.7|) then shows that a pole at ?/o = 
only appears for m = 2. However, there is only one diagram per order of this kind, which 
is shown in Fig. 9(b). Since this diagram factorizes into loops of the kind of Fig. 9(a), 
it contributes a pole in y^ of order — 1. Therefore the pole at order A^ = 2 is the 
only primitive singularity in the series (|A.3| ) for the free energy; analogous arguments 
apply to the expansions of the contact field correlation functions. It follows that the 
renormalization can be carried out to all orders; see also refs. []59|, BO] . The loops of Fig. 9 



form a geometric series. Summing this series defines the ^-factors |43 



1 Up , 

2/0 



Zr 



dwo 
duf 



1 Up 

yo J 



(A.ll) 



and the flow equation iip 
function of Mp, 



y^up — cMp. The Casimir amplitude ACp becomes a regular 



ACp(Mp) 



v{v - 1) 



Up + 0{u\ 



(A.12) 



The coupling up defined by ( |A.11D is not unique. Any family of diffeomorphisms up ^ 
{up,yo) with fixed point u'p{0,yo) = will preserve the regular functional dependence 

{c/Co)up 



(A.13) 



Uj 

of observables such as ( |A.12| ). In particular, the linear transformation up 
leads to the flow equation ( p.39| ) and the Z-factors 

7 ^ ^0 ^ 7 ^0 (, Co y 

Zp = — 1 Up , Zp = — 1 Up 

Co \ yo J \ yo ) 

We may also substitute Co(?/o) by Co(?/o = 0) in ( |A.13|) and (|2.39| ) so that the dependence 
of the flow equation on yo is contained only in the linear term. This scheme has been 
called "minimal subtraction" in ref. However, the Z-factors ( |A.13|) still have a 

regular part as yo ^ generated by the function ciyo). Eliminating this part requires 
a non-linear transformation of Mp and introduces cubic and higher order terms into the 



flow equation Of course, these terms are spurious; i.e., they do not generate higher 
order terms in expressions like (|2.40|) for local observables. (In the present case of a single 



coupling constant Mp, the value of the Co drops out of these expressions as well. In more 
general cases with several coupling constants, ratios of operator product coefficients do 
play a role; see the discussion in ref. 0].) 

The replica trick is unproblematic within perturbation theory, since it reduces to con- 
venient bookkeeping of the averaging over disorder. Eqs. ( A.13|) , ( p.39| ) are independent of 
p, and in Eq. (|A.12|) , the dependence on p reduces to the combinatoric prefactor. Hence, 
the random limit AC = Ymip^oCp/p is trivial and leads to Eqs. ( |3.62| ) and ( p.63| ). 

The exponents (|3.64| ) can also be obtained in a different way. They follow from the 
fact that the two-string interaction (jp^cj? does not renormalize the "mass term" (jx^ at any 
order. Consider the (normal-ordered) density field $(to) = 00(r = O,to)5 which has the 
two-point function 

°° ( — a ]^ 

mtoMt',));{uo,R)-^ ^ 



Af=0 



A^! 



dtoi . . . dtow($(to)$(to)$2(toi) • • • $2(t0^)): 



(A.14) 
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Its short-distance asymptotics is related to the return probabihty of strings to the origin 
r = 0. In the linear theory, 



mtoMQ); \to - tx'^^ (A.15) 

for I to ~ ^ol/-R^ ^ 1- Any perturbative correction to this exponent arises from the renor- 
malization of the fields $(to) and $(to). The renormalization of $(to) is due to a short- 
distance coupling of the form 

$(to)$2(toi) . . . ^2{toN) = t^'^'^^Clfin, . . . , rjv_i) $(to) + . . . (A.16) 

for to/R'^ (with toj - t = toTj, t^j - % = UTjk for j, A; = 1, . . . , A^, and n, . . . , tn_i 
denoting a basis of the fixed ratios Tj,Tjk), and there is a corresponding expression for 
$(to). However, it is obvious that the product on the l.h.s. couples only to contact fields 
of at least two lines, and therefore C^f = at all orders A^. The singularity ( A.15|) remains 



unchanged, C* = Co = 1/2- The strong-coupling form ( |3.86| ) of the singularity cannot be 
reached by perturbation theory. 

Now we turn to the dynamic path integral. After the change of variables to ~^ ^oto, 
ho {vq/ a'^Y^'^hQ, and Iiq {ol/voY^'^hQ, ( [4.3|) takes the normal form 



j VhoVhoexp - J dvdt I -hi + ho ldt,ho - h/^ho 




(A.17) 



with Aq = —go- As discussed in ref. |0, the transformation ([4.17|) leads back to ( |A.1| ) 



rendering dynamical perturbation theory identical to replica perturbation theory. The 
path integral ( |A.1|) has to be supplemented with boundary conditions corresponding to 
the random limit p 0. In this limit, the boundary condition (|A.2| ) describes the one- 
point-function of an initially flat surface, {h{r,to = 0)) = 0. At any integer value of p, the 
dynamical correlations are artificially truncated, (/i(ri,t) . . . h{rm,t)h{r, L))p = for all 
m > p by (lOI). 

In the context of growth dynamics, the representation (|A.17|) has an advantage over 



( |A.1| ): it contains explicitly the most interesting observable, the height field h. The price 



to pay is that perturbation theory becomes more complicated. The differences to replica 
perturbation theory are two-fold: (i) There is the cubic vertex hCVhy. (ii) The vertex 
h"^, the image of the replica interaction 0^0^, no longer has incoming lines. These vertices 
can therefore be integrated out, giving the standard diagrammatics with the two kinds of 
propagators ( ^I9| ) and ( [4.10| ). 



We now discuss the renormalization of ( |A.17| ) and show its equivalence to replica renor- 
malization up to one-loop order. The real-space response function has the diagrammatic 
expansion shown in Fig. 10(a). To order Uq, the expansion reads 

i?'^(/io(r,to)/io(r,to + i?'))(Mo) =i?'(/io(r,to)/io(r,to + i?'))(0)+ 0(i/o%o,M^) • (A.18) 
The one-loop diagram does not have a pole at = 2 since its short- distance singularity 



cancels with a geometric factor 2 — d [0, [T^, ^ . Constructing the strong-coupling fixed 



point in d = 1 requires taking into account the resulting finite renormalization of the 
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(a) 




(b) 




Fig. 10: Diagrammatic expansions generated by the dynamic functional ( f4.3D . The hnes with one 
and two arrows denote the unperturbed response function (4.9) and the unperturbed correlation 
function ( [4.10| ), respectively. Dots represent the vertices /i(V/i)^; each incoming line to a vertex 
has to be differentiated with respect to r. (a) Response function (/i(ri, toi)^(r2> ^02))- The one- 
loop diagram is regular at yo = 0. (b) Stationary growth rate (3tg/io)_R. The boxed subdiagram 
contains a simple pole at yo = 0. A further diagram with a loop as in (a) is regular at yo = 
and has been omitted. 



response function JT^. For the critical fixed point above d = 2, however, it can be 
ignored. 

The expansion for the growth velocity is shown in Fig. 10(b). The tadpole diagram at 
order Aq has the form 

Xo[I{R) - lim I{R)] = R-'^ (A.19) 

with 

J(i?) = I dr'dt', [VMr',t',)h{r,t))n{uo = 0)]' ; (A.20) 

the second term in ( |A.19| ) is generated by the normal-ordering ( [4.15| ) and cancels the 
ultraviolet divergence of I{R)- At order Aq, the boxed subdiagram contributes a pole 
originating from the integration region where the two vertices approach each other. Hence, 
we have to this order 

R'-'^'idtMRiuo) = -1 fl - —uo + 0{yyo,ul)] . (A.21) 

The pole can be absorbed into the definition of the variables hp = Zp^hQ, hp = Zp^ho, 
tp = Zptto, Up = ZpUo, with 

Zphiup) = 1 + ^up + Oiul), (A.22) 

Zptiup) = 1+0(4), (A.23) 

Zp(up) = l-—up + Oiul). (A.24) 

yo 

This reparametrization respects ( [4.25| ) and renders both the response function and the 
growth rate regular as ?/o ^ 0, 

R''{hp{r,tp)hp{r,tp + R'')){up) = R''{ho{r,to)ho{r,to + R'')){0) 
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+ 0{yQUp,ul) , 



(A.25) 



R'-'''{dt,hp)ji{up) = -l4/2(i + o(|/oV,4)) • 



(A.26) 



It leads to the beta function (|2.38| ) and to the exponents 

X* = 2 - z* = xo - Up log Zph{up) = 0{yl) 



(A.27) 
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